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Abstract 


This paper used the first-principle calculations based on the density functional theory 
(DET) to determine the electronic structure and optical absorption of rutile TiO; co-doped with 


nickel (Ni) as a transition metal and cerium (Ce) as a rare earth metal. 


The calculated results demonstrated that the nickel and cerium TiO: co-doping induced 
impurity energy levels above the valence band maximum (VBM). There reduced transition 
electrons required energy and limit the recombination rate of charge carriers, as well as 


prolonging their life span. 


Likewise, the optical spectra of Ti0: co-doped with nickel and cerium revealed a larger 


photo-response for visible light than that of nickel or cerium, single doped TiO:. 


Keywords: Ni-Ce, co-doping, rutile TiO;, DFT. 


Résumé 


Ce manuscrit constitue des calculs de premiers principes basés sur la théorie de la 
fonctionnelle de la densité (DFT) pour déterminer la structure électronique et l'absorption 
optique de l’oxyde de titane TiO>, dans sa phase rutile, co-dopé avec du nickel (Ni) comme métal 


de transition et du cérium (Ce) comme un métal des terres rares. 


Les résultats de calculs ont démontré que le co-dopage de TiO; au nickel ou au cérium 
induit des niveaux d'énergie d'impureté au dessus du maximum de la bande de valence (VBM). 
Ceci réduit l'énergie nécessaire aux électrons de transition et limite le taux de recombinaison des 


porteurs de charge, tout en prolongeant leur durée de vie. 


De même, les spectres optiques du TiO, co-dopé avec le nickel ou le cérium ont révélé 
une plus grande photo-réponse pour la lumière visible que celle de TiO: dopé seulement au nickel 


ou au cérium. 


Mots clés : Co-dopage Ni-Ce, TiO; rutile, DET. 
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Introduction générale 


Les oxydes semi-conducteurs ont reçu une attention considérable en tant que solution 
possible pour l'énergie propre et pour la dégradation photocatalytique des polluants dangereux 
pour l'environnement. Depuis 1970, de nombreux efforts ont été réalisés pour explorer ces 
matériaux et optimiser leurs configurations de bande d'énergie pour des applications 
potentielles. Parmi ceux-ci, citons l’oxyde de titane Ti; qui est utilisé dans les cellules photo 
électrochimiques car 1l a la capacité de photo-dissocier l'eau et de produire du photo-courant 


pour la production d'hydrogène. 


TiO; peut photo-dégrader de nombreux polluants organiques dissous dans l'eau par 
irradiation UV. Cependant, en raison de sa grande bande interdite, il absorbe mal la lumière 
visible et est donc inefficace pour la photocatalyse dans le domaine du visible. De nombreuses 
tentatives ont été faites pour améliorer l'activité photocatalytique de TiO, sous irradiation aux 
UV et / ou à la lumière visible. La stratégie principale consistait à diminuer l'énergie de la 
bande interdite, c'est-à-dire des énergies de l'ordre de 1,6-3,0 eV, en particulier à 2,5 eV, 


intensité maximale d'émission. 


Trois caractéristiques principales sont requises pour obtenir une bonne activité de la 
lumière visible pour TiO:: i) le dopage doit introduire des états dans la bande interdite 
pouvant être excités en absorbant la lumière visible, ii) le minimum de la bande de conduction 
(CBM) doit être au-dessus du niveau de réduction H: / H:,0 et le maximum de bande de 
valence (VBM) ainsi que les états dopants doivent être inférieur au niveau d'oxydation H20 / 
O; pour assurer une activité photo réductrice, 1ii) recouvrement suffisant entre les états 
d'énergie de dopant et les états de bandes de TiO; est nécessaire, afin de permettre un transfert 


rapide des porteurs photo-excités à la surface pendant leur durée de vie. 


TiO: fut alors l'objectif de beaucoup de recherches depuis quelques années jusqu'à 
maintenant. Généralement, le dopage est un moyen pour modifier et améliorer les propriétés 
électroniques et optiques de T10, [1-4]. Citons le dopage avec les non-métaux (N, C, B, S) [5- 


16], les métaux de transition [17], le co-dopage peut être aussi utilisé. 


Dans des travaux de Yin et al. [18] ont étudié une série d’éléments donneur-accepteur 
pour modifier les structures de bande du TiO: afin d'améliorer la dissociation photo 
électrochimique de l'eau. Ils ont trouvé que les non-métaux carbone et azote (C et N) sont les 
meilleurs candidats en tant qu'accepteurs, alors que les métaux de transition Mo, W, Nb et Ta 


sont les meilleurs candidats en tant que donneurs. Cependant, ceux ci en tant que dopants 


peuvent réduire la mobilité des porteurs par la formation d'états fortement localisés dans la 


bande interdite. 


D’autres recherches ont montré que deux dopants ont un effet plus synergique qu'un 
seul dopant, en ce qui concerne l'augmentation de l'absorption dans la région de la lumière 


visible et l'amélioration de l'activité photocatalytique du TiO;, [19-23]. 


Une meilleure efficacité photocatalytique peut être obtenue si Ti02 est co-dopés avec 
un métal de transition et un métal des terres rares. Ce travail consiste à co-dopé TiO: par : un 
métal de transition le nickel Ni et un métal des terres rares le cérium Ce pour mieux 
comprendre les différents phénomènes dus aux effets de dopage. La sélection de ces deux 
métaux est basée sur les résultats trouvés pour chacun d’eux pris comme un seul dopant 


[20,24-26]. 
Le manuscrit est réparti en trois chapitres organisés comme suit : 


Chapitre I. Il est constitué de trois parties principales. Dans la première partie, nous 
développerons une étude bibliographiques détaillée sur l’oxyde de titane, ses caractéristiques, 
ses propriétés et ses applications particulièrement en photocatalyse. Dans la deuxième partie 
on présentera les fondements de la théorie de la fonctionnelle de la densité (DFT) et l’accent 
sera porté sur la partie échange et corrélation de l’énergie, qui conditionne de façon cruciale la 
qualité des résultats obtenus. La troisième partie sera consacré à la description du formalisme 
de la méthode des ondes planes augmentées linéarisées à potentiel total ( ful potentiel 


liniarised augmented plane wave : FP-LAPW) utilisée ainsi que l’algorithme du code de 


calcul WIEN2K. 


Chapitre IL. Dans ce chapitre, nous étudions les propriétés structurales et électroniques des 


composés Ti1-x À x; O2 (x = 0, 0,25, 0,50 et 0,75). (A= Ce, Ni). 


Chapitre IIL Dans ce chapitre sont présentés les résultats des propriétés structurales, 


électroniques et optiques des composés T1o,50Ce025Ni0 2505. 


Enfin, les principaux résultats acquis de la présente étude sont résumés dans une conclusion 


générale. 
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Chapitre 1 


Revue bibliographique 


Ce chapitre est une étude bibliographique détaillée de l’oxyde de titane, ses caractéristiques, 
ses propriétés et ses applications ainsi que les fondements de la théorie de la fonctionnelle de 
la densité (DET) et le formalisme de la méthode des ondes planes augmentées et linéarisées à 


potentiel total (FP-LAPW). 


Chapitre 1 : Revue bibliographique 


Partie A : L’oxyde de titane dans la littérature 
L INTRODUCTION 
II. PROPRIETES STRUCTURALES 
III. PROPRIETES ELECTRONIQUES 
IV. LES PROPRIETES OPTIQUES 
V. DOPAGE DU DIOXYDE DE TITANE 
V.1. Dopage cationique 
V.2. Dopage anionique 
V.3. Co-dopage de TiO> 
VI. PHOTOCATALYSE HETEROGENE 
VL.1. Introduction 
VI.2. Principe 
VL.3. Photocatalyseurs aux TiO2 
VII. LES APPLICATIONS DE TiO; 


Partie B : Théorie de la fonctionnelle de la densité DFT 
L. INTRODUCTION 
II. EQUATION DE SCHRÔDINGER 
IL. APPROXIMATION DE BORN-OPPENHEIMER 
IV. APPROXIMATION DE HARTREE-FOCK 
V. APPROXIMATION DE THOMAS-FERMI 
VI. THEORIE DE LA FONCTIONNELLE DE LA DENSITE (DFT) 
VL.1. Théorèmes de Hohenberg et Kohn 
VI.2. Les équations de Kohn-Sham 
VL3. La fonctionnelle d’échange-corrélation 
VI.3.a. Approximation de la densité locale (LDA) 
VI.3.b. Approximation du gradient généralisé (GGA) 
VI.4. Résolution des équations de Kohn-Sham 


Partie C : Méthode des ondes planes augmentées et linéarisées FP-LAPW 
L. INTRODUCTION 
II. METHODE DES ONDES PLANES AUGMENTEES LINEARISEES (FP-LAPW) 
IL.1. Méthode des ondes planes augmentées (APW) 
IL.2. Principe de la méthode FP-LAPW 
IL.3. Les rôles des énergies de linéarisation E; 
IL.4. Développement de la LAPW en orbitales locales 
11.4.a. La méthode LAPW+LO 
11.4.b. La méthode APW+lo 
ILS. Le concept de la méthode FP-LAPW 
III. LE CODE WIEN2K 
IIL.1. Vue d'ensemble du code 
IIL.2. Détails de calcul 


Partie A : L’oxyde de titane dans la littérature 
L INTRODUCTION 


L’oxyde de titane ou TiO, est un matériau populaire qui a de nombreuses applications dans 
différents domaines [1]: il sert habituellement de substrat de support pour les particules 
métalliques en catalyse hétérogène; sa stabilité en fait un matériau de revêtement 
anticorrosion parfait; ses propriétés diélectriques en font un matériau important dans les 
appareils électroniques; la nature semi-conductrice en fait également un capteur de gaz 
potentiel; 1l est également utilisé comme un bon pigment blanc dans la peinture et les 


cosmétiques pour son indice de réfraction élevé. 


TiO, a été introduit dans le domaine de la photocatalyse depuis longtemps car il est 
responsable du "farinage" des peintures, le complexe organique dans la peinture peut être 
photo décomposé à partir de TiO:. Cependant, il n'est pas devenu l'étoile de la photocatalyse 
avant que Fujishima et Honda ne l'utilisent pour décomposer avec succès H0 en 1972 [2]. 
D'une part, d'énormes efforts ont été consacrés à améliorer son efficacité énergétique ; d'autre 
part, il a été largement appliqué dans les domaines du traitement des eaux usées et de la 
purification de l'air [3], de l'auto-nettoyage [4] et de l'anti-brouillard [5] où les processus 


photocatalytiques joueraient un rôle important. 


Dans cette première partie de ce chapitre, nous passons en revue les structures géométriques 
et les propriétés électroniques et optiques du TiO;, ainsi que l'effet du dopage particulièrement 
sur les propriétés optiques et le gap énergétique. Ceci est suivi des principes de base et des 


processus clés de la photocatalyse sur TiO». 


IL. PROPRIETES STRUCTURALES 


L'oxyde de titane TiO2 appartient à la famille des oxydes de métaux de transition. Il existe 
trois polymorphes connus dans la nature qui sont l'anatase (nom dérivé du mot grec anatasis 
signifiant extension), le rutile (nom dérivé du latin rutilus = rouge) et le brookite (nommé en 


l'honneur du minéralogiste H.J. Brooke) [6,7] 


Les trois polymorphes sont tous composés de chaînes d'octaèdres TiO4 de base reliés 
différemment par des arêtes et des coins. Dans l'anatase de structure cristalline tétragonale, 
chaque octaèdre partage des coins et quatre arêtes avec les quatre autres octaèdres, ce qui 
donne une structure en Zigzag (Fig. la). Pour le rutile, aussi de structure cristalline 


tétragonale, deux bords d'octaèdres sont partagés. Les chaînes de TiO, sont liées ensemble par 


8 


une liaison en coin pour former une structure linéaire (Fig. 1b). Dans le cas des brookites, les 
octaèdres partagent trois bords et l'arrangement octaédrique produit une structure cristalline 


orthorhombique (fig. 1c) [8]. 


Dans le rutile, chaque octaèdre est en contact avec 10 voisins (deux partageant une arête et 
huit partageant un coin), tandis que dans l’anatase, le nombre de coordination est 8 (quatre 
partageant une arête et quatre partageant un coin). L'anatase et le rutile sont tous deux 
tétragonaux avec respectivement 12 et 6 atomes par maille unitaire. Dans chacune des 
structures, chaque atome de titane est lié à six atomes d'oxygène et chaque atome d'oxygène 
est lié à trois atomes de titane. La maille unitaire du rutile contient deux atomes de Ti à (0,0,0) 


et (72, V2, 2) et quatre atomes d'oxygène à (u, u, 0) et (2 +u , 2 -u, 2) 


La symétrie octaédrique (Où) de l'octaèdre TiO6, dans le rutile, est déformée selon la symétrie 
D:r. cette distorsion est due à la présence des différentes longueurs de liaison Ti-O dans le 
plan (équatorial) et hors du plan (axial), ainsi que de la déviation des deux types d'angles de 
liaison Ti-O-Ti de 90°.Dans l’anatase, la distorsion de l'octaèdre est davantage due aux 
déplacements supplémentaires des atomes d'oxygène, générant une symétrie DA locale. Les 
distances de Ti-O sont plus courtes dans l’anatase (1,91 À et 1.95 À) que dans le rutile (1.94 
À et 1.99 À), mais les distances inter ioniques Ti-Ti sont plus grande dans l'anatase (3,79 À et 
3,04 À) que dans le rutile (3,57 À et 2,96 À). La structure cristalline de la brookite diffère de 
celle de l'anatase et du rutile par les longueurs de liaison Ti-O et les angles de liaison Ti-O-T1. 
où il existe six liaisons Ti-O différentes allant de 1,87 À à 2,04 À et 12 angles de liaison 
différents allant de 77 à 105°.La brookite est un polymorphe plus rare que l'anatase et le rutile 
et aussi plus difficile à produire en phase pure. Elle forme communément une phase 


secondaire mineure avec l'anatase et le rutile [9-13]. 


Le rutile est la forme thermodynamiquement stable dans la plupart des conditions de 
température et de pression et, de ce fait, est moins soluble dans l'HF et le H,SO, concentré 
que l'anatase ou la brookite [14]. Cela ressort également des calculs théoriques 
thermodynamiques. Cependant, il est indiqué que l'anatase est la phase stable lorsque la taille 
des particules est suffisamment petite. Ceci a été expliqué par une énergie de surface de 
l’anatase inférieure à celle du rutile lorsque la taille des particules est nanométrique[15]. Le 
recuit de l’oxyde de titane à diverses températures conduit aux transformations de phase 


suivantes (Fig.2) : anatase à brookite à rutile, anatase à rutile, brookite à anatase à rutile. 


Figure 1 : Structures cristallines de TiO; dans la phase : (a) rutile, (b ) anatase, (c) brookite. 
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Figure 2: Transformations de phase dans TiO: en fonction de la température de recuit [16] 
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Ces transformations de phase dépendent des méthodes de préparation, de la taille des 
particules et du dopageégalement. La transition de la phase oxyde de titane amorphe à la 
phase anatase se produit autour de 300°C. La transformation de la phase rutile se produit 
lorsque la température de recuit dépasse 700°C et également que la transformation du rutile 
est un processus irréversible. Hanaor et al [17] ont décrit la transformation de l'anatase en 
rutile en faisant varier la température de transition par rapport à la pression. La figure 3 
montre les transitions de phase dans TiO> par rapport à la température de recuit et à la 
pression de travail. Comme il est rapporté, la transformation de la phase anatase en rutile se 


produit autour de 605° € sous une pression de 1 atm. 


n° 


c 
x 


+ 
- 


O/Ti ratio 


Figure 3 : Diagramme de phase du dioxyde de titane adapté de la référence [18] 


Un autre facteur affectant la transformation de phase est la présence de défauts sur le sous- 
réseau d'oxygène (augmentation des lacunes d'oxygène TiO:,). Des lacunes d'oxygène sont 
également connues pour être créées dans l'atmosphère d'hydrogène, favorisant ainsi la 
transformation en rutile. Lorsque la phase rutile et la phase anatase possèdent toutes deux une 
structure tétragonale, le rutile est plus dense et possède donc une plus grande densité (4,26 g / 
cm°) que l'anatase (3,89 g / cm). Les films de TiO> amorphe possèdent une large gamme de 
densités, de 2,4 g / cm” pour les films poreux [19] à 3,2 - 3,65 g / cm°[20]. Les films de TiO; 


ayant une densité plus faible favorisent la diffusion des impuretés, donc dans le cas du dopage 
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avec des impuretés [21]. Jetons un coup d'œil sur la figure 3 qui montre les équilibres de 
phase Ti-O du dioxyde de titane par rapport au rapport oxygène/titane (O/T1) et à la 
température. La figure révèle que TiO2 est un oxyde pauvre en oxygène Ti0,.Dans le 
diagramme de phase, comme nous le voyons, il y a un grand nombre de composés qui varient 
en fonction de la stœchiométrie. A faible taux stœchiométrique, on peut trouver l'évolution de 
Ti, TO, TiO, Ti203, Ti3Os. A un rapport stœchiométrique supérieur entre Ti305 et TiO, il 
existe une évolution de la série homologue Ti:O,1, appelée phases de Magneli ou phase 


cristallographique (n > 4) [18]. 


Une étude approfondie a été adoptée dans la littérature basée sur l'effet des impuretés et des 
dopants sur la transformation de la phase anatase en rutile. Il a été rapporté que l'inclusion 
d'une certaine quantité d'impuretés ou l'addition de dopants modifient remarquablement la 
transformation de phase. Il a été rapporté aussi que ces dopants peuvent également agir 
comme inhibiteurs ou promoteurs de la transformation de phase en fonction de la taille, de la 


valence, etc. 
Le tableau 1 résume les paramètres physiques des différentes structures de l’oxyde de titane 


Tableau 1 : Caractéristiques des différentes structures de Ti0; [22] 


propriétés rutile Anatase brookite 
Taille moyenne de cristal 0.17 à 0.24 um 0.14 à 0.17 um / 
Couleur Bleu à crème bleu Brun foncé au noir 
verdâtre 
La masse volumique (g.cm°) 4.26 3.84 4.17 
La masse moléculaire relative 79.9 79,5 79.9 
dureté (échelles de mohs) 5.5à6.5 5.5 à 6.0 5.5 à 6.0 
La liaison Ti- O(À) 1.959 1.917 1.9 
La liaison Ti- Ti(À) 2.96 / / 
La structure quadratique quadratique Orthorhombique 
Les paramètres de maille (À) | a=b=4.549 | a =b = 3.784 c-9.514 a =9.184 
c = 2.959 b=5.447 
c= 5.145 
Densité volumique (Kg/m:) 4240 3830 4170 
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III. PROPRIETES ELECTRONIQUES 


La détermination de la structure des bandes électroniques, les densités d’états totales et 
partielle (TDOS, PDOS), le gap d’énergie joue un rôle important pour la connaissance des 


propriétés électroniques des matériaux. 


TiO> se présente comme un matériau semi-conducteur avec une large bande interdite 
légèrement supérieure à 3 eV. Dans son état stœchiométrique, les états occupés dérivent 
principalement des atomes d’oxygène, les états libres dérivent des atomes de Ti. Une fois 


réduit, le matériau contient des défauts telles que les lacunes d’oxygène. 


Les énergies de gaps du rutile, d’anatase et de la brookite ont pour valeur 3 eV, 3,2 eV et 3,1 
eV respectivement [23,24]. Ces valeurs de gap donnent lieu à des transitions correspondant à 


des photons se situant dans le domaine de l’ultraviolet. 


La structure électronique de l'anatase et du rutile est étudiée par différentes approches 
théoriques et expérimentales telles que la théorie fonctionnelle de densité, la théorie de la 
densité hybride, la spectroscopie de perte d'énergie électronique. . [25-28]. Les calculs de la 
densité des états (DOS) révèlent que le bas de la bande de valence O 25 est situé à 1,7 eV, 
tandis que le haut est O 2p avec une largeur de 5,17 eV [25,26]. Le bas de La bande de 
conduction (CB) est à 8 eV et se compose principalement d'états T1 3d, qui présentent deux 
structures distinctes, en dessous et au-dessus de 5 eV. Les bandes de conduction au-dessus de 


8 eV sont principalement de caractères s et p. 


Pour l'analyse des états de liaison deTiO; rutile, la partie supérieur de la bande de valence est 
composée d’orbitale O2p d’une largeur de 6,22 eV. La partie inférieure est la bande O2s ayant 
une largeur de 1,94 eV. L’intervalle énergétique entre O 25 et le minimum de CB est de 17,98 
eV. La partie inférieure de la bande CB est constitué de deux états Ti 3d d’une largeur de 5,9 
eV [291]. Le diagramme des orbitales moléculaire et le calcul électronique révèlent que Ti 3d 
est divisé en deux orbitales e, ett.. Les atomes d’oxygène entourent trois atomes de titane et, 
par conséquent, l’hybridation sp” des atomes d’oxygène forme trois liaisons 6 dans le plan et 
une liaison r perpendiculaire au plan. Deux orbitales Ti se combinent avec l'orbitale Ti sp” ét 
forment six liaisons © pointant vers les six atomes d'oxygène et forment l'orbitale moléculaire. 
Dans le cas de l'anatase, la DOS est aussi décomposé en Tics, T1 t28. Ops et Opr. Le haut de la 
bande de valence est divisé en trois régions. La région inférieure est la liaison 6 où Ops 


contribue principalement à la liaison, la liaison x dans la région de l’énergie moyenne et les 
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états non liés O p, dans la région supérieure. La bande de conduction est décomposée en 
Tie>5 eV et t<5eV. Le minimum de la bande de conduction est d,,. Le reste des bandes t2 
est lié aux états x. Il y a très peu de rapports sur la structure électronique du TiO> en phase 


brookite [30,31]. L'image des orbitales moléculaires de TiO, rutile et anatase est montrée sur 


la figure 4. 
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Figure 4 : Image des orbitales moléculaires de TiO; (a) anatase (b) rutile [32,33]. 
IV. PROPRIETES OPTIQUES 


Lorsque la lumière, d'une énergie suffisante, interagît avec un matériau, elle induit une 
transition électronique,des états occupés, situés sous le niveau de Fermi vers les états 


inoccupés situé au-dessus du niveau de Fermi. 


Il existe deux principaux types de transition de bande interdite dans les semi-conducteurs. 
Pour les semi-conducteurs à bande interdite directe, la transition électronique de la bande de 
valence à la bande de conduction est autorisée, l’absorption électronique et l’émission sont 
généralement fortes. Pour les semi-conducteurs à bande interdite indirecte, la transition 
électronique de la bande de valence à la bande de conduction est interdite et la transition est 
assistée par phonon, ce qui implique que l'énergie et la quantité de mouvement de la paire 


électron-trou sont modifiées dans la transition. 


TiO présente à la fois une bande interdite directe et indirecte pour l’anatase et le rutile 
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respectivement. Les propriétés optiques de ces phases dépendent de la disposition des atomes 
de titane et d’oxygène dans les octaèdres de Ti06 [34]. La transition interdite directe est due à 
la transition L'(VB) — L' (CB), alors que la transition indirecte autorisée est due à la transition 


M — F de O 2p à Ti 3d dans la zone de Brillouin. 


La bande interdite directe pour le rutile apparaît à 3,031 eV et celle indirecte à 3,049 eV. La 
transition indirecte dans le rutile est assistée par une énergie de phonon acoustique de 11 meV 
et une énergie de phonon optique de 15 meV. L’intervalle énergétique entre la bande interdite 
directe et indirecte est estimée à 10-50 meV. Comme dans la phase rutile, le gap indirect de 


l'anatase TiO; est dû à la transition M — T' de O 2p à Ti 3d. 


En régime nanocristallin, la bande interdite de l'anatase T10: est décalée vers le bleue par 
rapport à sa valeur de 3,2 eV. Cette augmentation de la bande interdite est principalement 


attribuée à l'effet de confinement quantique de la taille réduite du TiO:. 


Les propriétés optiques sont quantifiées et généralement décrites par la constante diélectrique 
€ (w), par l'indice de réfraction n (@), le coefficient d'extinction k (œ@) et le coefficient 


d'absorption à (@). 


TiO: est un matériau ayant un indice de réfraction élevé (tableau.2) et parmi les trois phases 
cristallines stables, c’est le rutile qui possède le plus grand indice. Ceci combiné à un 
coefficient de diffusion de la lumière visible élevé, font de la phase rutile un pigment blanc de 
choix pour l’industrie (peintures, colorant alimentaire ou pharmaceutique) [35,36]. Aussi ses 
propriétés optiques ainsi que sa non toxicité lui permettent d’être employé comme écran anti 


UV pour les crèmes solaires [37,38]. 


Tableau 2:quelques propriétés optiques de TiO: 


Phase Rutile Anatase Brookite 
Densité [39] 4,13 10 3,99 

Indice de Réfraction[40] 2,605-2,903 | 2,561-2,488 | 2,583-2,700 
Energie du gap /eV [41], [34] 3 3,20 3,14 
Module d’Young / GPa[421], [43] 243 194 199-245 
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V. DOPAGE DU DIOXYDE DE TITANE 


Une des applications les plus importantes du TiO, est la photocatalyse. Mais, en raison 
de la largeur de bande interdite et de l'absorption uniquement dans la région UV, ce semi- 
conducteur ne peut pas être utilisé dans la photocatalyse à la lumière visible. Le dopage est 
une méthode permettant de créer des états intermédiaires dans la bande interdite de TiO, afin 
de déplacer le bord d'absorption vers la zone visible d'éclairage. L'incorporation d'impuretés 
augmente la capacité d'absorption de la lumière visible en réglant la bande interdite de TiO; 


pur. 


La réponse optique d’un matériau est déterminée par sa structure électronique, elle- 
mêmeentièrement liée à la structure cristallographique, la composition chimique (du volume 
et de lasurface) ainsi qu’à la dimension du matériau en question (confinement pour les 
dimensionsnanométriques).Le dopage cationique (remplacement de l’ion Ti**) ou anionique 
(remplacement de l’ion O7)peut modifier les propriétés optiques du TiO2. Dans ce cas il est 
souhaitable de conserver lastructure cristallographique initiale pour préserver les propriétés 
photocatalytiques du TiO».Du fait du rayon ionique et de la charge de l’anion O”, il est plus 


facile de substituer le cationTi** avec d’autres métaux de transition. 


La structure de bande des semi-conducteurs peut être modifiée par le processus de dopage. 
Pour le TiO, le procédé est habituellement mis en œuvre afin d'étendre sa plage de réponse 
optique dans la région visible en remplaçant les ions Ti** ou O dans la structure par des 
substituant métalliques ou non métalliques. Les sites de dopage ne sont pas limités à la 
substitution d'atomes dans le réseau mais peuvent également être situés entre les sites Ti et 
O7,en positions interstitielles [44,45ST, Les espèces dopantes sont généralement incorporées 
dans la structure pour réduire la valeur Eg en abaïssant le bord supérieur de la bande de 
conduction ou en élevant le bord inférieur de la bande de valence. Cela a pour effet de 


sensibiliser le Ti0, aux longueurs d'onde de la gamme visible. 
V.1. Dopage cationique 


Différents métaux ont été employés pour doper le TiO: par des méthodes de chimie en 
voie humide [46], de traitements à haute température, ou d’implantation d’ions [47]. D’après 
lalittérature, on peut citer des dopages avec les ions métalliques V, Cr, Mn, Fe, Co, Ni, Cu, 
Zr,Sn, W,... les lanthanides La, Ce, Nd,... [48] ou encore avec les alcalins Li, Na, K [491]. 


Selonle pourcentage de cations substitués, plusieurs propriétés vont être modifiées : la photo- 
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réactivité et le taux de recombinaison des charges [46], la stabilité thermique, la tailledes 
cristaux, la cristallinité, ou encore la concentration de défauts de surface [50]. La 
structureélectronique s’en trouve donc modifiée via la densité d’états. Des calculs de bandes 
basés surla théorie de la fonctionnelle de la densité (DFT) ont montré la formation d’un 
niveauélectronique occupé dans la bande interdite du TiO: dont la position varie selon le 
dopant et sa concentration [51]. D’autres auteurs ont montré la réduction du gap du TiO: de 
0,55 eVgrâce à la substitution avec des ions Nd** et la création d’un nouveau niveau non- 
occupé audessus de la BV du TiO, [52], ainsi que le comportement de semi-métal de TiO; 
dopé Nd ou V grâce à la présence des états électroniques Nd 4f et V 3d en bas de la BC 
[53,54]. 


Chen et al. [55], à l'aide de la spectroscopie d'absorption X, ont examiné la réduction 
de la bande interdite de TiO, lors du dopage du cérium. Il a attribué le rétrécissement de la 
bande interdite à la présence de défauts oxygénés et d'états 4f associés à Ce** / Ce** situés 
sous la bande de conduction. Ces bandes entraînent une réduction efficace de la bande 
interdite. Zhao et al. [56] ont démontré théoriquement la réduction de la bande interdite de 
TiO2 dopé avec différents types d'éléments de terres rares. Ils ont rapporté que les états 
électroniques 4f ne subissaient pas d'hybridation avec le niveau O2p ou Ti 3d du TiO:. Au 
lieu de cela, ces états 4f forment des états électroniques peu profonds près du bord de la bande 
de conduction, ce qui favorise la transition électronique de la bande de valence à ces états 4f. 
Ils ont calculé que lors du dopage du cérium et du néodyme, la bande interdite de TiO; 


diminuait respectivement de 0,12 eV et 0,07 eV. 


Expérimentalement, les procédures de synthèsepeuvent avoir un effet sur la structure de bande 
résultante. Elles appartiennent généralement à l'une des trois catégories suivantes : chimie en 
solution, traitement à haute température ou implantation ionique. À titre d'exemple, Anpo et 
al. Ont rapporté une différence significative dans les propriétés optiques résultantes du Ti0O; 
dopé au chrome (Cr) préparé par une méthode d'imprégnation humide comparée à celles 
synthétisées par implantation ionique [57]. L'augmentation de la réponse optique s'est avérée 
plus élevée pour le matériau synthétisé par la méthode d'implantation ionique où un décalage 
de bande vers la bande rouge a été observé. Dans la même publication, le groupe a utilisé la 
procédure d'implantation d'ions pour étudier la réponse optique du TiO: dopé avec une 
gamme de métaux de transition. Ils ont trouvé un décalage vers le rouge correspondant aux 


différents dopants tel que V> Cr> Mn> Fe> Ni. 


LE 


Un autre facteur clé dans l'utilisation des dopants dans les photocatalyseurs à base de TiO, est 
leur capacité à jouer le rôle de pièges à électrons et à trous, retardant ainsi la vitesse de 
recombinaison des porteurs de charge. Choi et al. ont publié une étude complète des effets du 
dopant et ont aperçu que la caractéristique la plus importante était les propriétés de transport 
de charge améliorées du TiO; résultant d'un dopant efficace. [58]. Le groupe a proposé que 
les dopants agissent à la fois comme pièges à électrons et à trous afin d'être photo-actifs. En 
effet, l'électron ou le trou piégé peut se recombiner avec son homologue mobile. Cependant, 
des dopants avec des configurations électroniques telles que Li*, Mg”, Al et Zn”*, se sont 


avérés avoir une influence insignifiante sur la photo-activité de ces systèmes. 


Finalement, ce sont les propriétés optiques qui vont être bouleversées avec, en général, 
un décalage vers le visible dans le seuil d’absorption [46], [48], [51], [52]. Encore une fois, 
lavaleur du décalage peut varier avec la concentration de dopant. Ce décalage du gap est 


attribué à une transition électronique entre le niveau électronique d du métal et la BC du TiO:. 


En général, le dopage de TiO: avec des ions métalliques dépend du type et de la concentration 
du dopant ainsi que du mode opératoire de son introduction et doit être considéré avec soin en 


fonction de l'application du photocatalyseur. 
V.2. Dopage anionique 


La plupart des éléments non métalliques tels que B, C, N et F ont également été utilisés 
comme dopants de TiO2:.Une série de techniques comprenant la synthèse sol-gel, le dépôt par 
pulvérisation cathodique et le dépôt chimique en phase vapeur, avec des études approfondies 
sur l'activité de la lumière visible d'une large gamme de composés à base de T10: dopés avec 


des élémentsnon métalliques sont réalisés [59,601]. 


Parmi les dopants anioniques, non métalliques, N s'est peut-être révélé le plus prometteur en 
termes d'amélioration de la photo-activité de TiO; dans le visible. [61] L'azote peut être 
facilement introduit dans le réseau de TiO, car il est similaire en taille atomique à l'oxygène, 
est très stable et a une faible énergie d'ionisation [62].Les premières études du dopage par 
l'azote a été effectué par Sato en 1986, lorsqu'il a rapporté la réponse à la lumière visible d'un 
photocatalyseur au TiO;, préparé en utilisant du NH4OH pendant l'étape de calcination. 
Cependant, ce n'est que lorsque les travaux d'Asahi et al. en 2001, la recherche sur le dopage 
du TiO; avec des ions non métalliques s'est intensifiée. Le groupe de travail a proposé un 


ensemble d'exigences pour obtenir une activité de lumière visible pour le TiO> à partir 
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d'études théoriques de DOS, soutenues par des observations expérimentales[63].Ils ont 
indiqué que, pour obtenir une activité de lumière visible de TiO», les critères suivants 
devraient être respectés: le dopant devrait produire des états dans le TiO02Eg absorbant la 
lumière visible; les nouveaux états d'impureté et le nouveau rebord CB inférieur résultant 
devraient être aussi élevés que ceux du TiO); les nouveaux états de l'Eg devraient se 
chevaucher avec les états de bande de TiO: pour permettre le transfert de porteurs de charge 


photo générés vers les sites actifs à la surface [63]. 


Comme avec les photocatalyseurs au TiO: dopés au métal, le procédé de préparation, et le 
type et la concentration des dopants non métalliques influencent la structure de la bande 


résultante et l'activité de la lumière visible du catalyseur. 
V.3. Co-dopage de TiO> 


Les propriétés de T10: peuvent également être modifiées par co-dopage, ce qui peut 
être obtenu par la combinaison de paires métal / métal, non métal / non métal ou non métal / 
métal. Un effet synergique (ie le non-métal peut provoquer un décalage rouge dans la région 
de la lumière visible, alors que le métal peut faciliter le transfert des porteurs de charge photo 
générés supprimant ainsi la recombinaison) peut être obtenu avec une combinaison appropriée 
de co-dopants TiO, dopé ou non dopé aux ions. Mis à part les co-dopants facilitant le 
rétrécissement des bandes interdites, leur combinaison peut entraîner la formation de 
différentes hétéro structures (c'est-à-dire différentes structures électroniques) par rapport au 
TiO:. Une hétéro structure constituée de différentes combinaisons de non-métal et de métal a 
la capacité d'améliorer la séparation des charges (métal) et l'absorption de la lumière visible 


(non-métal). 


Plusieurs combinaisons métalliques et non métalliques telles que N-I [57], C-V et Ag- 


V [54] ont été utilisés dans la photo dégradation de polluants et de colorants. 


Le co-dopage des ions lanthanides avec d'autres éléments est de plus en plus populaire en 
raison de leurs différentes propriétés d'absorption et de leur activité photocatalytique 
améliorée par rapport au T10: dopé individuellement [64]. Les métaux des terres rares ont été 
dopés avec d'autres ions de métaux des terres rares, des ions métalliques et des non-métaux, 
pour améliorer la réponse de la lumière visible à l'activité photocatalytique [65]. Le co-dopage 
des lanthanides avec des non-métaux tels que N, C, B sont largement étudié ces dernières 


années [66,67], tel que Ce / B-Ti0: [68] ,Ce / C-TiO: [69] et Er / Pr-T10: [70]. 
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Les résultats de ces chercheurs ont montré l'activité photocatalytique accrue et l'absorbance de 
la lumière visible par les combinaisons de métaux par rapport aux systèmes de TiO2 non 


dopés et dopé individuellement. 
VI PHOTOCATALYSE HETEROGENE 
VL.1. Introduction 


Plusieurs semi-conducteurs ont été testés dans la photocatalyse comme TiO, ZnO, 
ZnS, WO3, GaP, Fe203 et CdS. Les semi-conducteurs CdS et GaP ont l’avantage d’absorber 
par rapport au TiO; une fraction plus importante du spectre solaire, mais malheureusement 
ces semi-conducteurs sont instables durant la réaction. TiO), au contraire, est stable et 


constitue un bon compromis entre efficacité et stabilité. 
VI.2. Principe 


La photocatalyse hétérogène est un processus qui consiste à irradier, naturellement par le 
soleil ou artificiellement par lampe, un semi-conducteur. Ceci conduit à l’accélération de la 
photo réaction en faisant intervenir des réactions entre les paires électron/trou et les produits 
adsorbés à la surface du semi-conducteur. Lorsque ce dernier (SC) absorbe des photons 
d’énergie supérieure à la valeur de sa bande interdite (hv>Eg), il subit une excitation par 
rayonnement qui permet à un électron de la bande de valence d’être éjecté dans la bande de 


conduction, créant un site d’oxydation (trou h+) et un site de réduction (électron e-) suivant la 


réaction suivante : 
Semi — Conducteur — egc + hiy (1) 
OHias F hèv — OHias (2) 


VL.3. Photocatalyseurs aux TiO> 


Les premiers effets photo actifs du TiO: ont été rapportés par Lenz en 1920 qui a observé un 
assombrissement du matériau dû à l'irradiation UV [71]. Goodeve et Kitchener ont ensuite 
décrit ses effets photo sensibilisants en 1938 lorsqu'ils l’ont utilisé comme pigment blanc dans 
les peintures et les papiers peints [72]. Cependant, ce n'est qu'en 1972 avec le travail de 
Fujishima et Honda que les recherches intensives sur les propriétés photocatalytiques de TiO2 


ont commencé sérieusement. 


Au cours des dernières décennies, les photocatalyseurs au TiO, ont été utilisés dans une large 
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gamme d'applications, de la décomposition de l'eau et de la conversion du combustible 
solaire, à la décomposition des polluants organiques pour l'assainissement de l'environnement 


[73-77]. dans les matériaux pour surfaces antimicrobiennes et auto-nettoyantes [78-80]. 


ES 


TiO, est un semi-conducteur qui passe à un état d'énergie élevée en recevant de l'énergie 
lumineuse, et libère des électrons de sa surface éclairée. Si l'énergie reçue à ce stade est 
suffisamment élevée, les électrons initialement situés dans la bande de valence sautent tous 
dans la bande de conduction. La paire de trous d'électrons est générée, ce qui génère des 
radicaux libres capables de subir une réaction secondaire. Cela forme le principe de l'activité 
photo-catalytique. Le trou produit a un fort pouvoir oxydant et les électrons ont un fort 


pouvoir réducteur. 


Le processus photocatalytique repose sur l’excitation du TiO2 par un rayonnement 
lumineux de longueur d’onde inférieure à 400 nm. Un électron passe de la bande de valence à 


la bande de conduction: 
TiO, -— TiO, + héy + ec (3) 


Les trous h: réagissent avec des donneurs d’électrons tels que l’eau, les anions OH adsorbés et les 
produits (R) adsorbés à la surface du semi-conducteur en formant des radicaux hydroxyles et des produits 


intermédiaires notés R° selon les équations suivantes : 


H; Ouas DEA H* + OH£as (4) 
OHoas + h? a OHas (5) 
Rads + he HA Robe (6) 


Les électrons réagissent avec des accepteurs d’électrons tels que le dioxygène pour former des radicaux 


super oxydes. Cette réaction limite la recombinaison des charges. 
O;,+e 0; (7) 


En l’absence d’accepteur et de donneur d’électrons appropriés, on assiste à l’annihilation trou/électron c’est 


la réaction de recombinaison très rapide de l’ordre de la picoseconde: 


TiO; + häy + epc — TiO) (8) 
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La figure 5 résume le mécanisme photocatalytique à la surface de TiO; 


Figure 5 : Schéma du mécanisme de dégradation photocatalytique à la surface de TiO: [811]. 
L'activité photo-catalytique d'un semi-conducteur est donc largement contrôlée par : 


a. Les propriétés d'absorption de la lumière comme le spectre d'absorption de la lumière et le 


coefficient d'absorption. 
b. Taux de réduction et d'oxydation à la surface par électrons et trous. 
c. Taux de recombinaison des trous d'électrons. 


d. Une surface importante avec une densité de surface absorbante constante conduit à une 


vitesse de réaction photo-catalytique plus rapide. 
VIL APPLICATIONS DE TiO; 


Comme tous les oxydes de métaux de transitions, le dioxyde de titane présente de 
nombreuses applications pouvant être classées en quatre catégories relatives à ses propriétés 
optiques, photocatalytiques, électriques et ses propriétés d’hydrophobicité. Etant facilement 
réductible, le dioxyde de titane occupe une place privilégiée parmi ces oxydes. Parmi les 
applications, son utilisations dans les systèmes photocatalytiques, les phénomènes "super- 
hydrophile", les cellules solaires à colorant, les guides d’ondes et les systèmes électro 


chromes. 


Le TiO, est également un opacifiant efficace sous forme de poudre, où il est utilisé 


comme pigment pour apporter de la blancheur et de l'opacité à des produits tels que peintures, 
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revêtements, plastiques, papiers, encres, aliments, médicaments (pilules et comprimés) ainsi 


que la plupart des dentifrices. 


Dans les produits cosmétiques et de soin, le dioxyde de titane est utilisé comme 
pigment, crème solaire et épaississant en raison de son indice de réfraction élevé, de sa forte 


capacité d'absorption UV et de sa résistance à la décoloration sous ultraviolets. 


En particulier sous la forme d'anatase, c’est un photo-catalyseur sous la lumière 


ultraviolette (UV) le plus efficace et le plus respectueux de l'environnement [82]. 
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Partie B : Théorie de la fonctionnelle de la densité DFT 
I. INTRODUCTION 


La théorie de la fonctionnelle de la densité (DFT, ou Density Functional Theory) est l'une des 
méthodes, largement utilisées pour les calculs ab initio de la structure des atomes, des 


molécules, des cristaux, des surfaces et de leurs interactions [1, 2]. 


L'idée principale de la DET est de décrire un système interagissant à plusieurs corps via sa 
densité de particules et non via sa fonction d'onde. Sa signification est de réduire les 3N 
degrés de liberté du système à N-corps à un système qui a seulement trois coordonnées 
spatiales à travers sa densité de particules. Elle est utilisée pour calculer les propriétés d’un 
solide telles que: les propriétés électroniques, optiques, les géométries d'équilibre, les 


énergies d'activation, les énergies de réaction. . 


Dans cette partie du chapitre, nous discuterons de la théorie générale et du contexte historique 
de la DFT. Nous abordons ses aspects de ses implémentations pratiques ainsi que quelques 


concepts importants. 


II. EQUATION DE SCHRÔDINGER 


L'équation de Schrüdinger [3] est le point de départ de toute discussion de la mécanique 


quantique. La forme indépendante de l'heure de cette équation est 
ÎY =EY (1) 


Où À est l'opérateur Hamiltonien, E est l'énergie totale du système et W est la fonction d'onde 


du corps. 


Considérons un système formé par M noyaux et N électrons. Les particules telles que les 
électrons et les noyaux interagissent les unes avec les autres, c'est-à-dire que le potentiel de 
coulomb et, par conséquent, l'équation de Schrôdinger en unités atomiques peuvent s'écrire 


comme suit : 


N 1 M 1 N M y NN 1 M M D 7 
.-ÿin-Pn-DD MN N EE € 
2. * 2M4 i c A AT; RAB @) 


Où RA et ri sont les positions de M (noyaux) et N (électrons), respectivement. La distance 


entre le noyau A et l'électron i est ri = | ri - RA |. La distance entre les électrons i et j est : 
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rj ={|ri-r;|et la distance entre les noyaux A et B est Ras = | Ra - Re |. ZA est le numéro 
atomique du noyau, tandis que MA est sa masse en unités de la masse d'un électron. V4 et 


V£ sont les opérateurs laplaciens. 


Dans l'équation ci-dessus, les deux premiers termes représentent respectivement l'énergie 
cinétique des électrons et des noyaux. Le troisième terme est l'attraction de Coulomb entre les 
noyaux et les électrons, tandis que la répulsion entre les électrons et les noyaux est 


représentée respectivement par les quatrième et cinquième termes. 


IL. APPROXIMATION DE BORN-OPPENHEIMER 


L'équation de Schrôdinger peut être résolue facilement pour l'espèce moléculaire la plus 
simple (H*?), mais elle est très compliquée à résoudre dans le cas des systèmes moléculaires 
plus complexes. Dans ce cas, l'approximation de Born-Oppenheimer [4] est utilisée, elle joue 


un rôle important dans les calculs de structures électroniques. 


ES 


Cette approximation consiste à séparer l'équation de Schrôdinger en parties nucléaires et 
électroniques. Elle est basée sur le fait que la masse d’un noyau quelconque est nettement 
supérieure à celle de l'électron. Ainsi, dans la plupart des cas, les noyaux peuvent être 
considérés comme des points fixes et les électrons se déplacent sous le potentiel externe fixe 
produit par les noyaux. La fonction d'onde électronique dépend donc uniquement de la 
position des noyaux et non de leur mouvement. Sous cette approximation, le second terme de 
l'équation (2) est négligeable ainsi que la répulsion entre les noyaux, le dernier terme peut être 
considéré comme une constante pour les noyaux fixes. Les termes restants sont le 


Hamiltonien électronique (Éarec) qui peut être écrit comme 


N N M NN 
Ban) 3) dt dd ® 
“ 26 a Ar; 
La fonction d'onde électronique totale pour la molécule est donc : 
Fatec Wélec = Eelec Velec (4) 
et la solution générale de l'Eq (2) peut être finalement écrite comme : 
Yi (noyaux, électrons) = Ÿ (électrons) Ÿ (noyaux). 


L'énergie totale est évidemment la somme de l'énergie électronique et nucléaire. L'énergie 
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électronique est la somme de l’énergie potentielle et de l'énergie cinétique des électrons. Les 
électrons se déplacent dans le champ électrostatique généré par les noyaux et sont soumis à la 


répulsion électronique: 


AZB 
Etot = Eetec + s » R (5) 
4 AB 


=1 B>A 


Dans ce cas, la résolution de l'équation (2) est encore complexe pour la fonction d'onde à 


électrons multiples, contenant 3N variables. 


Afin de résoudre ce problème de plusieurs électrons, des approches spéciales ont été 


développées, telles que la théorie de Thomas-Fermi et la théorie de Hartree Fock. 
IV. APPROXIMATION DE HARTREE-FOCK 


La méthode de Hartree-Fock a été introduite à la fin des années 1920 et développée en 1927 
par D. KR. Hartree. C'est une approximation pour la détermination de la fonction d'onde et de 
l'énergie d'un système à plusieurs corps dans un état stationnaire [5]. Cette théorie simplifie 
l'interaction électron-électron, en élargissant la fonction de nombreuses ondes électroniques 
en un produit de la fonction d'onde électronique unique. La fonction d'onde de l'état 
fondamental d'un système à N électrons est écrite comme un seul déterminant de Slater à N 
dimensions, Wsp, avec les fonctions d'une onde électronique (Wy;) comme composantes. Le 


déterminant de Slater est : 


; WiCX1) Wi1CX2) ++ Wi1(Xn) 
Vel] = [00 HO TG © 
Un(X1) Un (X2) Le Un(Xn) 


Avec cette fonction d'onde, l'équation complète de Hartree-Fock peut être écrite comme : 


- LUE V., +®, kr (r)+ Va P(r)= (7) (7) 


Où Vexen Vi (r) est le terme non-local d’échange ajouté 
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Vu P(r)=-X far EE Y(r) (8) 


fl 
ik) Ir r| 


Elle se compose de quatre termes. Les deux premiers termes représentent la contribution de 
l'énergie cinétique et potentielle des électrons. Le troisième terme est également connu sous le 
nom de terme de Hartree qui est le potentiel électrostatique généré par la distribution de 
charge des N électrons. Le quatrième terme, ou terme d'échange, résulte de l'exclusion du 


principe de Pauli. 


L'approximation Hartree-Fock suppose une image à un seul électron de la structure 
électronique, c'est-à-dire que la distribution des N électrons est donnée par la somme des 
distributions à un électron |yl. Cependant, la théorie de Hartree-Fock a été établie en 
supposant une forme à un seul déterminant pour la fonction d'onde et elle néglige également 
la corrélation entre les électrons. Par conséquent, le potentiel moyen non local des électrons 
provenant des autres électrons n'est pas pris en compte, ce qui conduit à une mauvaise 
description de la structure électronique. Bien que la théorie HE soit qualitativement correcte 
pour divers matériaux et composés, elle ne permet pas de faire des prédictions quantitatives 


précises. 


La fonction d'onde électronique d’HF est construite en utilisant des fonctions d'onde 
électronique antisymétriques. Puisque dans les méthodes HE, les corrélations électroniques 
sont négligées et produisent avec de larges bandes passantes, des énergies plus élevées et des 
largeurs de bande trop faibles, l'approche HF est particulièrement adaptée aux petites 
molécules organiques et aux oxydes, qui ont un petit nombre d'électrons localisés. Cependant, 
cette approximation HF est moins appropriée pour les métaux de transition à densité 


électronique élevée. 


En raison de l'interaction de la densité de charge moyenne des électrons, l'énergie de l'état 
fondamental calculée dans la théorie HF est normalement plus élevée que l'énergie réelle de 
l'état fondamental. La théorie échoue dans les systèmes d'électrons délocalisés parce que cette 
méthode ignore le criblage collectif de Coulomb. Suite à cette faiblesse de la théorie HP, la 


méthode DFT a été développée, qui traite de la densité électronique. 
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V. APPROXIMATION DE THOMAS-FERMI 


La théorie de Thomas [6] et Fermi [7] est une véritable théorie de la fonctionnelle de la 
densité (DFT) en ce sens que toutes les contributions à l’énergie électronique totale (cinétique 


aussi bien qu’électrostatique) sont exprimées en termes de la densité électronique. 


ES 


Elle consiste à subdiviser le système inhomogène en volume élémentaires dr (soit 
implicitement dv), dans lesquelles les électrons ont un comportement d’un gaz homogène de 


densité constante. Afin de rester dans cette limite, le potentiel V4 doit être une fonction de r 


. \ 1957 , . 1 , à : 4 
variant lentement à l’échelle de la longueur d’onde de Fermi mo C est-à-dire devant resté dans 
F 


le domaine des vecteurs d’onde inférieurs à la norme de kr. 


+ SU Re a er 
Le potentiel est défini par l’action des noyaux (V,, (F) —— — ) ainsi que par l’action 


électrostatique des électrons en r2 sur ceux en r;. 


S'agissant 1c1 de deux électrons 1 et 2 : 


Veff — Vext (r) + Vee ( l1,12) (9) 


La complexité de la solution de Thomas-Fermi provient du terme de répulsion inter- 


électronique et le terme d’énergie potentielle dû à V.,. est donné par : 
E, = [o(rh(rd?r (10) 


En effet, Vee (r1, r2) donne un terme bi-électronique difficile à calculer : 
V. = jj 2) d°r, (11) 
Ir = ;| 


Pour ce faire le remplacement de ce terme par une intégrale de Coulomb traduisant 
l’interaction de deux densités électroniques a été tenté par Thomas et Fermi. Cette intégrale 


s’écrit : 
1 
Zaol= fe) ptdr, (12) 
12 


C’est un pas de plus par rapport au terme direct de Hartree mais ça ne fut pas une issue 
satisfaisante pour la solution du problème car J72 /p] ignore l’influence réciproque de la 


position instantanée des électrons et donc les effets d’échange-corrélation. 
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VI THEORIE DE LA FONCTIONNELLE DE LA DENSITE (DFT) 


La DFT est la méthode la plus importante, elle est largement utilisée en raison de sa bonne 
précision et de son efficacité de calcul élevée. Celle-ci prend en compte la structure des 
solides cristallins et de leurs interactions. Cette méthode fonctionne sur la base de la 
distribution de densité de charge électronique plutôt que de la fonction d'onde à plusieurs 
électrons. L'avantage de l'utilisation de la densité électronique est une dimension beaucoup 
plus réduite par rapport à celle de la fonction d'onde. En effet, cette méthode est facile à 


appliquer aux systèmes plus importants car elle fonctionne avec des densités électroniques. 


Des discussions détaillées sur la DFT peuvent être trouvées dans de nombreux livres [8, 91 et 
dans des articles de revue [10-12]. Bien que la méthode DFT soit précise pour les états 
fondamentaux, elle ne tient pas suffisamment compte de la forme de la fonctionnelle 
d'échange et de corrélation. Certaines approximations sont donc utilisées pour affiner le 


modèle pour les interactions échange-corrélations. 
VL.1. Théorèmes de Hohenberg et Kohn 


Hohenberg et Kohn [13] ont présenté en 1964 la théorie moderne de la DFT des systèmes à 
plusieurs corps où la théorie fonctionnelle de la densité (DFT) est basée sur deux théorèmes 
simples, mais extrêmement importants, connus sous le nom de théorèmes de Hohenberg- 
Kohn. Nous présentons ici ces deux théorèmes comme indiqué dans le livre Electronic 


Structure : Basic Theory and Practical Methods, publié par Martin [14]. 


+ Théorème I. Pour tout système de particules en interaction dans un potentiel externe V,(r), 
le potentiel V,,{r) est déterminé uniquement, à l'exception d'une constante, par la densité de 


particules à l'état fondamental no(r). 


Puisque le potentiel et par conséquent le Hamiltonien sont déterminés, la fonction d'onde à 
plusieurs corps peut être déterminée pour l'état fondamental ainsi que pour tous les états 
excités. Par conséquent, il est alors conclu de ce théorème que toutes les propriétés du 
système sont complètement déterminées étant donné seulement la densité de l'état 


fondamental nor) [15]. 


+ Théorème 2. Une fonction universelle pour l'énergie E[n] en termes de densité n(r) peut 
être définie, valable pour tout potentiel externe V.,{r). Pour tout V,.{r) particulier, l'énergie 


exacte à l'état fondamental du système est la valeur minimale globale de cette fonction, et la 


Le 


densité n(r) qui minimise la fonctionnelle est la densité exacte de l'état fondamental no(r). 


De ce théorème, il résulte que les fonctions Ef[n] permettent de déterminer l'énergie et la 


densité de l'état fondamental [15]. En fait, cette fonctionnalité peut être considérée comme 
Eln] = F[n] + Î Vu (n(r)dèr + En (13) 


Où En est l'énergie d'interaction des noyaux, et F[n] est une fonctionnelle universelle de la 
densité incluant toutes les énergies internes (cinétique, T{n] et potentiel, E;. [n]) du système 


d'électrons en interaction [15]: 
Fin] = Tin] + Enn] (14) 


Alors que les théorèmes de Hohenberg-Kohn prouvent l'existence d'une solution exacte à 
l'état fondamental du système, 1ls n'offrent pas une approche pratique pour son calcul. 
Cependant, en 1965, environ un an après les travaux de Hohenberg et Kohn, ont proposé une 


méthode réalisable pour les calculs DET qui est décrite dans la section suivante. 
VI.2. Les équations de Kohn-Sham 


Pour appliquer le formalisme DET, il faut une bonne approximation de la fonctionnelle 
Ffn(r)]. L'approximation correcte de Kohm-Sham (KS) [15] consiste à remplacer le problème 
du système original à plusieurs corps par un système de particules indépendantes auxiliaires 
qui aura la même densité d'état fondamental mais sera plus facile à résoudre. L’avantage 
principal de cette approximation est qu'elle permet de calculer l'énergie cinétique de la 
manière la plus simple. Elle fournit clairement une image des systèmes électroniques 
interagissant avec une particule. Cette méthode a été introduite en 1965 et Kohn a reçu le prix 


Nobel de chimie en 1998, avec John Pople. 


L'Hamiltonien auxiliaire indépendant est : 


A 


1 
Hxs = _ + Vxs(r) (15) 


Afin d'évaluer l'énergie cinétique d'un système à N particules indépendantes, 1ls ont obtenu le 
potentiel correspondant VKs (r), en résolvant le système à un électron de N équations de 


Schrôdinger: 
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1 
(-57 Fe Hs) bi@r) = air) (16) 


Chacune des N orbitales wi (r) a un électron avec les valeurs propres les plus basses &i, avec la 


densité : 
N 
pG) = D IWC (17) 
i=1 
La condition de normalisation est 
N = [ ecer (18) 


L'énergie cinétique des particules indépendantes Ts [p(r)] des particules non-interactives est 


TLp(r)] = > [vi ovacar (19) 


La fonctionnelle universelle F[p(r)] est 


Flp(r)] = Tsle(r)] + Elo) + Exclp(n)] (20) 


Où Eg[p (r)] représente l'énergie électrostatique classique des électrons, 


GE De PE dr dr' (21) 


et Exc[p(r)] est l'énergie d'échange-corrélation. 


Vks(r) = VexeCr) + Var) + Vrc(r) (22) 


Où Va(r) est le potentiel de Hartree Vy(r) 


M = [ar (23) 


Et V,c(r) est le potentiel d'échange-corrélation donné par : 


6Excle(r)] 
Vxc(r) = en (24) 
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Les équations (22), (23) et (24) sont les équations de KS. Ces équations sont auto-cohérentes 


car Vxc (r) dépend de la densité. 


La définition de l'énergie fonctionnelle E,. est 


Exec) = Tle)] - Tslp()] + Einrlo@)] + Exlp(r)] (25) 


où T{p(r)] et Eix [p(r)] représentent respectivement les énergies cinétiques et d'interaction 


électron-électron du système. 
VL3. La fonctionnelle d’échange-corrélation 


Les équations ci-dessus fournissent une bonne méthode pour trouver l'énergie de l'état 
fondamental d'un système en interaction lorsque la forme de Exc est connue. Cependant, la 
forme de Exc a été calculée uniquement pour quelques systèmes simples, donc ce terme est 


généralement inconnu. 


Le terme Exc est couramment approché dans les calculs de structure électronique dans 
l'approximation de densité locale ou l'approximation de gradient généralisée. Au cours des 
dernières décennies, des centaines d'auteurs ont publié des formes différentes pour la fonction 
Exc. La première méthode a été proposée par Perdew et s'appelle aussi l'échelle de Jacob [16]. 


Il a présenté les types simples de fonction d'échange-corrélation qui sont les suivants : 
VL3.a. Approximation de la densité locale (LDA) 


La méthode la plus largement utilisée et la plus simple pour obtenir l'énergie d'échange et de 


corrélation s'appelle l'approximation de la densité locale (LDA : local density approximation). 


L’approximation de la densité locale LDA fait l’hypothèse que l’énergie d’échange et de 
corrélation est la même dans un solide que pour un gaz homogène d’électrons libres. L’idée 
de l’approximation LDA est de considérer le potentiel d’échange et de corrélation comme une 
quantité locale, définie en un point r, dépendant faiblement des variations de la densité autour 
de ce point. Ainsi, en un point r auquel correspond une densité p(r), il sera associé un 
potentiel d’échange et de corrélation comparable à celui d’un gaz d’électrons de même densité 


p(r). L'énergie d’échange - corrélation totale E,. s’exprime alors en fonction de l’énergie 


d’échange et corrélation par électron £,. telle que : 
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ERA [p(r)] = Î p(r)exc (otr))dr (26) 


En d’autres termes, on postule qu’autour de chaque point r, on peut remplacer le système réel 


par un gaz homogène d’électrons de densité p(r). La fonction £X°" (p ) peut être décomposée 


en une contribution d’échange £;(p) et de corrélation £.(p) : 


&(ptr)= 2 (etr)}+ £c(p(r) (27) 


Le terme d’échange d’un tel gaz a été déterminé exactement par la technique Monte-Carlo par 


Ceperley et al [17] qui ont tabulé le terme d’échange-corrélation ehom(b) en fonction du rayon 


de Wigner-Seitz : 


l 
3 
= Hi) (28) 
47p 
Il existe de nombreux travaux de paramétrisation de En ) comme par exemple ceux de 


Perdew et Zunger [18]. 


Il existe également une version de la LDA qui permet de prendre en compte le spin 
électronique : c’est l’approximation de la densité de spin locale LSDA. L’énergie d’échange 


et de corrélation Exc devient une fonctionnelle des deux densités de spin haut et bas : 


ERA op d)= [ons e Tr) ptr) br (29) 
avec 
Exc — ex” (o(r) (30) 


Cette approche fut initialement proposée par Slater et permet de résoudre certains problèmes 
liés à une approche LDA, notamment le traitement de systèmes soumis à des champs 


magnétiques et les systèmes où les effets relativistes deviennent importants. 


Le terme &xc est évalué à partir de calculs précis (Monte Carlo quantique) et par des 


procédures de paramétrisation. 


S 


Bien qu’étant une approche assez simple conceptuellement, l’approximation LDA permet 
néanmoins d’obtenir de bons résultats. La LDA à travers sa description locale de la densité 
donne des résultats d’une précision inattendue à première vue. Cependant certaines quantités 
telle que l’énergie de cohésion sont fortement surestimées, d’autres telles que les paramètres 
de maille et longueur de liaison sont sous-estimées. La LDA donne, par exemple une 
mauvaise estimation du gap des semi-conducteurs (40-50 %), ce qui n’est pas très surprenant, 
car cette grandeur ne relève que partiellement de l’état fondamental. Elle conduit parfois à une 
description médiocre des propriétés de certains systèmes ( l’état fondamental de Fe à la 
température ambiante est désigné comme étant la variété cfc y , alors qu’il s’agit en fait de la 


variété à cc ), surtout pour les actinides et les lanthanides. 


Cependant, le but poursuivi est la description des expériences et la reproduction des résultats 
obtenus. Cette théorie est un outil tant qualitatif (propriétés des liaisons chimiques) que 
quantitatif (moments magnétiques, champs hyperfins, ordres magnétiques...) assez 
performant, dont les lacunes avérées sont autant d’ouvertures sur différents champs 


d’investigation. 
VI.3.b. Approximation du gradient généralisé (GGA) 


L'objectif principal de l'utilisation de l'approximation (GGA: généralised gradient 
approximation) est que plusieurs systèmes sont souvent très différents d'un gaz d'électrons 
uniforme. Généralement, on peut dire que les systèmes réels (systèmes atomiques ou 
moléculaires) ne sont pas homogènes. Pour cette raison, le formalisme de correction pour 
décrire la fonctionnelle d’échange-corrélation a été amélioré. La GGA est une méthode non 


locale qui dépend uniquement de la densité électronique, E££A[p(#)] et son gradient V[p(r)] 


La fonctionnelle est de la forme [19]: 


Ejé" Lo(P)] = FAT on (32) 


Dans ce nouveau formalisme, le terme de l'énergie d'échange-corrélation par particule dans un 


système électronique inhomogène est représenté par : 


exe” [Le] VII] (33) 


Compte tenu de l'utilisation extensive de la DFT en tant que mécanisme pour décrire 
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différentes propriétés physiques des matériaux, différentes formules GGA ont été améliorées, 
l'approximation GGA a montré la fiabilité et la précision de la simulation pour différents 
matériaux. La fonction GGA améliore également la fiabilité de l'énergie totale [20,21] et les 
bandes d'énergie à la structure de la bande. Cette approche permet une meilleure description 


des systèmes magnétiques et optiques des systèmes non uniformes. 
La GGA proposé par Pedrew-Burke et Ernzerhoft [22] est utilisée dans cette thèse. 
VI.4. Résolution des équations de Kohn-Sham 


En utilisant des méthodes de particules indépendantes, les équations KS fournissent un moyen 
d'obtenir la densité et l'énergie exactes de l'état fondamental d'un système de matière 
condensée. Les équations KS doivent être résolues de manière cohérente car le potentiel KS 
effectif VKs et la densité électronique n(r) sont étroitement liés. Ceci est généralement 
effectué numériquement à l'aide de certaines itérations auto-cohérentes, comme illustré à la 
Fig.5. Le processus commence par une densité électronique initiale, généralement une 
superposition de densité électronique atomique, puis le potentiel KS efficace VKks est calculé 
et l'équation KS est résolue avec des valeurs propres et des fonctions d'onde, une nouvelle 
densité électronique est ensuite calculée à partir des fonctions d'onde. Après cela, les 


conditions auto-cohérentes sont vérifiées. 


Les conditions auto-cohérentes peuvent être le changement de l'énergie totale ou de la densité 
électronique de l'itération précédente ou la force totale agissant sur les atomes est inférieure à 
une petite quantité choisie ou une combinaison de ces conditions individuelles. Si la 
cohérence n'est pas atteinte, la densité électronique calculée sera mélangée avec la densité 
électronique des itérations précédentes pour obtenir une nouvelle densité électronique. Une 
nouvelle itération commencera avec la nouvelle densité électronique. Ce processus continue 
jusqu'à ce que l'auto-cohérence soit atteinte. Une fois l'auto-consistance atteinte, différentes 
quantités peuvent être calculées, dont l'énergie totale, les forces, les contraintes, les valeurs 


propres, la densité électronique des états, la structure de la bande, etc... 


L'étape la plus consommatrice de temps dans l'ensemble du processus consiste à résoudre les 


équations KS avec un potentiel KS donné par VKs. 


39 


Mélange 
de 
Pout (r) 


et 


Pin(r) 


Densité initiale 


Pin (r) 


Calcul du potentiel effectif de Kohen et Sham 


Vofr le] = Vext + Vale] + Vrc[p] 


Résolution des équations de Kohen et Sham 


EE + vs] ,(r) = e;(r)9,(r) 


Détermination de l'énergie de 
Fermi Er 


Calcul de la nouvelle densité électronique 


Poue(r) = D 1b:GIIE 


Calcul de 
l'énergie 


totale 
Converge _—_— - 


Figure 5 : Diagramme de la théorie de la fonctionnelle de la densité DFT. 
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Partie C : Méthode des ondes planes augmentées et linéarisées FP-LAPW 


L INTRODUCTION 


La théorie de la fonctionnelle de la densité est une approche puissante pour le 
traitement du problème à plusieurs corps. Cependant, 1l est important de faire le choix 
convenable d’une base de fonctions d’onde pour la résolution des équations de Khon-Sham. Il 
existe plusieurs méthodes de calculs des structures de bandes, leur point commun est la 
résolution des trois équations de Kohn et Sham de façon auto-cohérente. Ces méthodes 
sont classées en trois principaux types selon qu’ils nécessitent des résultats 


expérimentaux ou des données fondamentales : 
*Les méthodes empiriques pour lesquelles les calculs nécessitent des résultats expérimentaux. 


*Les méthodes semi-empiriques pour lesquelles les calculs nécessitant à la fois des résultats 


expérimentaux et des données fondamentales. 


*Les méthodes ab-initio pour lesquelles les calculs nécessitent seulement les données 


fondamentales. 


Ces dernières années, les chercheurs ont développé des méthodes ab-initio basées sur 
des concepts théoriques, c'est-à-dire sans paramètres ajustés sur l’expérience. Ces méthodes 
sont dites de premiers principes car leurs calculs sont faits à partir des premières informations 
qu'on a sur les interactions entre les électrons et les noyaux des atomes qui 


constituent le solide. 


Parmi ces méthodes on trouve la LCAO :Méthodes basées sur une combinaison 
linéaire d’orbitales atomiques [1,2], OPW : Les méthodes dérivées des ondes planes 
orthogonalisées [3,4], APW: Les méthodes cellulaires du type ondes planes 
augmentées[S]...etc. Les méthodes qui ont permis de gagner beaucoup de temps de 
calcul, sont les méthodes linéarisées de Andersen [6]: Ondes planes augmentées et 


linéarisées (LAPW) et orbitales « Muffin-Tin » linéarisées (LMTO). 
II. METHODE DES ONDES PLANES AUGMENTEES LINEARISEES (FP-LAPW) 


La méthode des ondes planes augmentées et linéarisées LAPW (linearized augmented 


plane wave), développée par Andersen [4], est fondamentalement une amélioration de 
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la méthode dite : méthode des ondes planes augmentées (APW) élaborée par Slater [6-91]. 
Pour bien comprendre la méthode LAPW, il faut savoir d’abord les principes sur lesquels 


repose cette méthode qu’elle en dérive. 
IL.1. Méthode des ondes planes augmentées (APW) 


La méthode APW est la technique la plus populaire pour la résolution de la structure 
électronique en utilisant les équations de Kohn-Sham (KS). Slater a proposé la technique en 


1937 [4, 6]. 


Dans cette approche, dans le schéma APW, la cellule unitaire est divisée en deux types 
de régions : (1) des sphères centrées autour de tous les sites atomiques r, de rayon R,, et (ii) la 
région interstitielle restante (Fig. 6). Dans ce cas, les fonctions d'onde sont étendues en onde 
plane (PW) dont chacune est augmentée par des solutions atomiques sous la forme d'ondes 
partielles, c'est-à-dire une fonction radiale multipliée par des harmoniques sphériques. En 
particulier, des solutions radiales de l'équation de Schrôdinger sont employées à l'intérieur de 
sphères non atomisées centrées sur les atomes et d'ondes planes dans la zone interstitielle 


restante. 


L'introduction d'un tel ensemble de base est due au fait que près des noyaux, les 
fonctions d'onde et le potentiel sont très similaires à celles d'un atome, tandis que les atomes 


sont plus loin. Les APW comprennent : 


Ainubn(T,E), T<Ry 


o()= ST sn. (35) 
ON? Ceexp(i(k + G).r), r el 
G 


est la fonction d'onde, Q est le volume de la cellule unitaire, r” est la sphère de position @ 
avec les coordonnées polaires r”, k” est un vecteur d'onde dans la zone de Brillouin 
irréductible (1BZ) et um est la solution numérique de l’équation de Schrôdinger radiale à 


l'énergie &. Les orbitales KS w (r) sont exprimées comme une combinaison linéaire de APW 


p(r). 


A l'intérieur de la sphère MT, une orbitale KS ne peut être décrite avec précision que 
si &, dans les fonctions de base APW, est égale à l'énergie propre ei. Par conséquent, un 


ensemble dépendant de l'énergie différent des fonctions de base APW doit être trouvé pour 
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chaque énergie propre. Cc et am Sont des coefficients de dilatation, Er est un paramètre (défini 
égal à l'énergie de la bande) et V la composante sphérique du potentiel dans la sphère. 


d? (+1) 
me HE + V(r) mn E; Tütm(T) = 0 (36) 


Slater a motivé l'utilisation de ces fonctions en notant que les ondes planes sont les 
solutions de l'équation de Schrôdinger dans un potentiel constant et que les fonctions radiales 
sont des solutions dans un potentiel sphérique. Cette approximation du potentiel est appelée 
"muffin tin" (MT). Puisque la continuité sur les limites des sphères doit être garantie sur la 
représentation double définie dans l'équation. (36), la contrainte doit être imposée. Dans la 
méthode APW, ceci est fait en définissant le um en termes de CG dans l'expansion harmonique 


sphérique des ondes planes. 


Région 
Interstitielle 


Figure 6: Représentation de potentiel Muffin-Tin «MT» 


Dans la technique APW, on continue à utiliser des ensembles de base (fonctions) qui 
couvrent tout l'espace dans la sphère. Cependant, APW est une technique couramment utilisée 
pour les calculs des propriétés structurales, électroniques, optiques et magnétiques des solides, 
mais même dans ce cas, elle présente certains inconvénients. Un exemple de la limitation de 
cette technique est: il ne peut pas être prolongé au-delà de l'approximation moyenne sphérique 
muffin-tin. Les fonctions de base, pour cette approche, ont un problème à la frontière du 
muffin-tin et donc à la frontière leurs dérivés sont discontinus. Un autre inconvénient de cette 
approche est que la fonction radiale um (r, E1) dépend de l'énergie, ce qui conduit à un 
problème de valeur propre non linéaire. Cela peut entraîner des complications numériques si 


uw devient très petit à la limite de la sphère vide. 
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IL.2. Principe de la méthode FP-LAPW 


La technique APW est modifiée par Anderson et est appelée méthode LAPW [5]. 
L'approche LAPW résout les problèmes de la méthode APW, c'est-à-dire que les fonctions de 
base et leur première dérivée sont interrompues à la frontière muffin-tin entre la région 
centrale et la région interstitielle. Pour résoudre le problème, Anderson a introduit un schéma 
de linéarisation et c'est pourquoi, il a utilisé le terme méthode LAPW (Linéarized Augmented 


Plan Wave). Il a rajouté un terme supplémentaire dans la base de la sphère MT. 


Dans la technique LAPW à l'intérieur de la région MT, la fonction d'onde radiale est 


linéarisée par une combinaison linéaire de Um (Tr, Er) et Wim (r, En [7-9] : 


D. au (re) + DE ü (r) Vu P) r<R, 
g(r) = 4 7" . (37) 
NES Cexpli(é + 6).F), r El 
G 


Où les bim sont des coefficients pour la dérivée d'énergie analogue à am. Les fonctions de 
base à l'intérieur des sphères sont des combinaisons linéaires de fonctions radiales and 1 (r), 
uw (r) et leurs dérivées d'énergie 1 (Tr) Yim (r) et wi (Tr) Yim (r) sont les fonctions 
d'augmentation. Les w sont définis comme dans la méthode APW selon l'équation (37) et la 
dérivée d'énergie, ui (Tr) Yim (r), vérifie l'équation suivante : 


d? IG+1 
D D + v6r) — Ei] rüm() = Tüm(r) (38) 


Les LAPW fournissent une base suffisamment flexible pour décrire correctement les 
fonctions propres avec des énergies propres proches de l'énergie de linéarisation, qui peuvent 
être maintenues fixes. Ce schéma nous permet d'obtenir toutes les énergies propres avec une 
seule diagonalisation contrairement à APW. Elles sont des ondes planes dans la zone 
interstitielle de la cellule unitaire, qui correspondent aux fonctions radiales numériques à 


l'intérieur des sphères avec l'exigence que les fonctions de base et leurs dérivées soient 


continués à la limite. 


Dans cette méthode, aucune approximation de forme n'est faite et par conséquent une 
telle procédure est souvent appelée "LAPW à potentiel total" (FP-LAPW). L'approximation 


beaucoup plus ancienne de muffin-tin correspond à ne retenir que les composantes 1 = 0 et m 
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= 0 dans l'équation (37). Une moyenne sphérique à l'intérieur des sphères et la moyenne 


volumique dans la région interstitielle sont ainsi prises. 


A l'intérieur de la sphère atomique est utilisée une combinaison linéaire de fonctions 
radiales multipliées par des harmoniques sphériques, Yim (r). La combinaison linéaire de w 
(r) et u(r) constitue la linéarisation de la fonction radiale. w (r) et u(r) sont obtenus par 
intégration numérique de l'équation de Schrôdinger radiale sur un maillage radial à l'intérieur 
de la sphère. Les LAPW ont plus de liberté variationnelle à l'intérieur des sphères que les 
APW. Cette plus grande flexibilité est due à la présence de deux fonctions radiales au lieu 
d’une. Les potentiels non sphériques à l'intérieur des sphères peuvent maintenant être traités 
sans difficulté. Il y a cependant un prix à payer pour la flexibilité supplémentaire des LAPW : 
les fonctions de base doivent avoir des dérivées continues et par conséquent des coupures 
d'ondes planes plus élevées un niveau de convergence donné. La solution des équations KS 


est élargie dans cette base combinée selon la méthode de variation linéaire : 
Pr = > Cn Prn (39) 
n 


Et les coefficients Cx sont déterminés par le principe variationnelle de Rayleigh-Ritz. La 
convergence de cet ensemble de base est contrôlée par un paramètre de coupure Rx X Kimax, 
où Rx est le plus petit rayon de sphère atomique dans la cellule unitaire et Kana est 


l'amplitude du plus grand vecteur K, dans l’équation (39). 
IL.3. Les rôles des énergies de linéarisation E; 


Les fonctions d’onde augmentées w et wi doivent satisfaire à la condition qu’elles 
soient orthogonales aux états de cœur dans la sphère MT. Mais cette condition n’est satisfaite 
que dans le cas où il n’y a pas d’états de cœur avec le même I. En conséquence, on risque de 
confondre les états de semi-cœur avec les états de valence. Le problème de non orthogonalité 
de quelques états de cœur n’est pas traité par la méthode APW, et l’utilisation de la méthode 
FP-LAPW nous oblige de choisir l’énergie E, délicatement. Alors, pour effectuer le calcul, E; 


doit subir une modification. 


La solution idéale dans de tels cas est d’utiliser un développement en orbitales locales. 
Cependant, cette option n’est pas disponible dans tous les programmes, et dans ce cas, on doit 


choisir un rayon de la sphère le plus grand possible. 
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Finalement, 1l faut remarquer que les divers E, devraient être définis indépendamment 
les uns des autres. Les bandes d’énergie ont des orbitales différentes. Pour un calcul précis 
de la structure électronique, E; doit être choisi le plus proche possible de l’énergie de la bande 


si la bande a le même 1. 
IL.4. Développement de la LAPW en orbitales locales 


Le but de la méthode LAPW est d’obtenir des énergies de bande précises au voisinage 
des énergies de linéarisation E;i [10]. Pour la plupart des matériaux, 1l suffit de choisir ces 
énergies au voisinage du centre des bandes. Ce qui n’est pas toujours possible car il existe des 
matériaux pour lesquels le choix d’une seule valeur de E; n’est pas suffisant pour calculer 
toutes les bandes d’énergie, c’est le cas pour les matériaux ayant des orbitales 4f [11,12] et les 
métaux de transition [13,14]. C’est le problème fondamental de l’état de semi-cœur qui est 
intermédiaire entre l’état de valence et celui de cœur. Pour pouvoir résoudre ce problème, on a 
recours soit à l’usage des fenêtres d’énergies multiples, soit à l’utilisation d’un développement 


en orbitales locales. 
I1.4.a. La méthode LAPW+LO 


Comme mentionné ci-dessus, dans certains cas 1l peut y avoir plus d'un état de valence 
(ou semi cœur) avec la même valeur de 1, ce qui rend difficile de choisir une énergie de 
linéarisation pertinente E1. Pour résoudre ce problème et améliorer la linéarisation, Singh et 
all.[15] ont introduit des orbitales locales (LO) qui augmentent la base LAPW définie pour 
certaines valeurs de 1. Cette méthode est connue sous le nom de LAPW + LO, où les LO sont 


définis comme 


T > Rur 


0 
— : 40 
pe) (latte Era) + BimÜi(r, Eva) + CmUi(r, Eu2) im) T < Rur (40) 


Ici E1 et E2 sont des énergies de linéarisation convenant respectivement aux états le plus 
élevé et le plus bas [16]. En d'autres termes, LO, o (r) est construit à partir des fonctions 
radiales LAPW à l'énergie E,1, et une troisième fonction radiale U à une seconde énergie E1 2, 
choisie pour améliorer au mieux la linéarisation. Puisque les LO ne sont pas connectés à des 
ondes planes dans la région interstitielle, ils n'ont pas de dépendance k ou G. Les trois 
coefficients Aïm, Bim €t Cim Sont déterminés par les exigences selon lesquelles les LO 


devraient être normalisés et devraient également avoir une valeur nulle et une pente à la limite 
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de la sphère muffin-tin [17, 18]. Par exemple, ajouter les LO pour les états p et d de tous les 
atomes d'une cellule avec N atomes augmente la base définie par (3 + 5) x N atome = 8N 
atome, ce qui n'est pas un grand nombre par rapport à la taille d’un ensemble de base typique 
[17]. Les LO se sont révélés plus efficaces pour améliorer la linéarisation que les méthodes 


alternatives avec des APW ayant des dérivées continues secondes et troisièmes [16]. 
11.4.b. La méthode APW+lo 


Les fonctions de base APW peuvent également être modifiées de la même manière que 
les fonctions LAPW en introduisant des orbitales locales. Sjosted et al. [191 ont amélioré 
l'APW en introduisant une base locale orbitale (lo), c'est-à-dire APW + lo. Comme dans 
l'approche APW, la fonction d'onde radiale est évaluée à énergie fixe mais cette nouvelle 
technique inclut un autre type d'orbitale pour la flexibilité. Ces orbitales sont notées lo au lieu 


de LO dans le LAPW [20]: 


6 = {le C0 + DOME 7 € Re 4 

La condition de normalisation que la fonction @ (r) aux rayons MT est nulle ne 
détermine les deux coefficients de l'équation. Dans cette approche, wi ne dépend pas des ondes 
planes et est inclus uniquement pour l'ensemble sélectionné des 1 nombres quantiques. Le 
terme dérivé d'énergie, dans la méthode APW + lo, est seulement inclus dans quelques lo et 
pas dans tous les plans vagues comme ça dans le LAPW. Par conséquent, nous ne sommes 
pas sûrs que la linéarisation d'énergie d’APW + bas soit précise comme LAPW, bien qu'elle 
converge plus rapidement que LAPW et donne la même précision que la technique LAPW 


[21]. 
ILS. Le concept de la méthode FP-LAPW 


L'approximation muffin-tin a été fréquemment utilisée en 1970 dans les techniques 
APW et LAPW [22]. Le potentiel est supposé constant dans la région interstitielle alors qu'il 
est sphérique symétrique dans la région muffin-tin. Ces approximations sont efficaces dans 


des systèmes hautement coordonnés comme les structures métalliques centrées sur le visage. 


Ces approches ne sont pas très efficaces dans les calculs des propriétés structurales et 
électroniques des solides à structure ouverte et à liaison covalente. Pour ce type de composés, 


les résultats calculés montrent une incohérence avec les résultats expérimentaux. Afin 
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d'obtenir de meilleures prédictions pour les propriétés de ses composés, aucune approximation 
de forme n'est utilisée. Ces composés sont traités avec la méthode FP-LAPW. Celle ci 
combine la base LAPW avec le traitement du potentiel total sans inclure l'approximation de 
forme dans les régions interstitielle et muffin-tin [23]. Dans l'approche FP-LAPW, le potentiel 
est étendu dans deux régions de telle sorte que (a) à l'intérieur de la sphère atomique :1l est 


dans les harmoniques de réseau et (b) dans la région interstitielle comme série de Fourier [23]: 


VG) = D Vau() Yu) (42) 
LM 

VOIS 2% exp(ik.r) (43) 
D LM 


Où (42) exprime le potentiel à l'intérieur de la sphère atomique, alors que (43) exprime le 
potentiel à l'extérieur de la sphère atomique. Cette forme est connue sous le nom de schéma 
de potentiel total. La sélection des rayons de la sphère dans l'approximation muffin-tin est très 


critique mais pas dans le FP-LAPW. 
III. LE CODE WIEN2K 
IIL.1. Vue d'ensemble du code 


Le Wien est un code qui fonctionne sous Linux ou UNIX, et 11 a été développé à Inst.f 
Materials Chemistry-Austria, publié par Blaha et al en 1990. Afin de l'améliorer et de le 
rendre adaptable à différents calculs et de traiter différentes propriétés physiques des 
matériaux, le code de Wien a subi plusieurs changements et chaque changement prend une 
notation publiée qui dépend de l'année de publication, par exemple Wien93, Wien95, Wien97, 
Wien2Kk (2000 année de publication). Cette thèse a été développée et réalisée avec la version 
Wien2k_14.2 qui présente le meilleur avantage en termes de temps de calcul et de formalisme 


utilisé pour le calcul des propriétés physiques des systèmes, l'efficacité et la fiabilité. 


W2web signifie le Wien to Web Le code Wien2k a démontré son efficacité dans le 
domaine de la chimie quantique et de la physique de la matière condensée, basé sur la 
méthode des ondes planes augmentées, mis en œuvre au sein de la DFT. Ce paquet est un 
ensemble de programmes indépendants écrits avec Fortran et calcule plusieurs propriétés 
physiques telles que: Propriétés électroniques des matériaux pour l'étude de la structure de 


bande, de la densité électronique totale / partielle (TDOS et PDOS) et des propriétés optiques. 
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Ce code permet de calculer les paramètres optiques tels que l'absorption, la réflectivité, 
l'indice de réfraction, la conductivité optique … etc. Il peut également traiter l'énergie totale du 
système et l'optimisation des structures, propriétés structurales, thermodynamique 
(enthalpie … etc), et propriétés magnétiques (calcul de l'état ferromagnétique et anti- 
ferromagnétique pour étudier la stabilité magnétique des systèmes), polarisation, gradients de 


champ électrique et champs hyperfins … etc. 
III2. Détails de calcul 


Nous commençons dans n'importe quel calcul des différentes propriétés physiques des 


matériaux par la procédure suivante : 


Premièrement, nous définissons le fichier “input” ou les données d'entrée avec 
l'extension case.struct, ce fichier montre la structure du matériau. Il contient les paramètres de 
la cellule unitaire, les positions atomiques du matériau et le nombre d'atomes différents. 
Ensuite, nous procédons à l'initialisation de notre calcul en exécutant une série de 
programmes qui seront entrés dans les programmes principaux du code, ainsi nous exécutons 


la commande init_lapw pour initialiser les sous-programmes suivants: 


NN : C'est un sous-programme qui initialise le calcul et détermine la distance entre les voisins 
les plus proches d'une certaine limite. Il a défini le facteur de distance F qui est souvent égal à 
2. Il traite et contrôle un chevauchement entre des sphères atomiques en ajustant le rayon 


atomique, son fichier de sortie est case.outputnn. 


Sgroup: Il définit et détermine le groupe d'espace de la structure étudiée, le fichier de sortie 


est case.struct.sgroup. 


SYMMETRY : C'est un programme qui traite les opérations de symétrie du groupe spatial 
étudié à partir des données d'entrée. Il détermine également le groupe de points de différents 
sites atomiques et la matrice de rotation correspondant aux opérations. Il est donc responsable 


de la génération des harmoniques sphériques Im. 


LSTART : Il est responsable de la génération des densités électroniques des atomes libres, et 
il nécessite la sélection du potentiel d'échange et de corrélation à utiliser dans le calcul : 
LSDA (Perdew et wang), GGA (Perdew et Wang et al Perdew 96 et 92). Dans ce sous- 
programme, nous spécifions l'énergie de coupure (coupure) qui sépare les états de valence et 


ceux du noyau, elle est généralement égale a -6,0 Ry. 
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Il est également à noter que dans le START, il peut montrer comment les orbitales sont 
traitées dans le calcul de la structure de bande, le traitement des états de valence avec ceux de 


semi-core ou core, et le calcul avec orbitale locale ou sans LO. 


KGEN : Il spécifie le nombre de K-point afin de créer une cellule dans la zone Brillouin 


uréductible (1BZ). 


DSTART : Il est basé sur une superposition de densités atomiques réalisées par le sous- 
programme LSTART, DSTART porte donc la densité initiale pour commencer à calculer un 


cycle auto-cohérent (SCF). 


Après avoir initié tous ces sous-programmes, tous les fichiers "input" requis dans le 
calcul SCF sont créés, et la SCF exécute un ensemble de programmes qui initialise la 
commande run_lapw, ce processus est itéré jusqu'à ce que la convergence (énergie, densité) 


soit obtenue. 


LAPWO0 : à partir de la densité, ce sous-programme calcule le potentiel de Coulomb [24] et le 


potentiel de corrélation d'échange, il génère le fichier case.scf. 


LAPW1 : calculer les valeurs propres et les vecteurs d'énergie en diagonalisant la matrice, 


ceci permet de calculer les bandes de valence. Le fichier de sortie est case.scfl et case.vector. 


LAPW2 : à partir des vecteurs propres et de l'énergie de Fermi, la densité de valence est 


calculée. Le fichier de sortie est case.scf2. 
LCORE : Il génère des énergies et des états de core et le fichier de sortie est case.scfc. 


MIXER : dans ce sous-programme, les densités des états de cœur, des états de semi-noyau et 
de valence sont ajoutées afin de produire une nouvelle densité totale qui sera une densité 


initiale pour une itération suivante. 


Ci-dessous nous présentons le diagramme de fonctionnement du code WIEN2XK : 
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SYMMETRY 


NN LSTART 
Vérifier le non Calcul atomique Fichier structure, [ DSTART 
Mon Hu = Enr Pyn fichier d’entrée 


Superposition des 
densités atomiques 


Ÿ 
Densité atomique Lomme 
fichier d’entrée — _KGEN 
Génération de la Loc 
maille k | 
LAPWO 


V?V. = —8np Poisson 
Vxc(P) LDA 


V=V + Vr 


Vur 


LCORE 


2 = 
E V2 + VIWx = Ex Pr Calcul atomique 


| Ex [ w Hu = Er Pr 
| Pcore [ Ecore 


| Pold 


MIXER 
Pnew — Pold @ (Pyai n.P ER | 


OUI _ NON 
STOP +——"""" CONVERGE 


Figure 7 : L’organigramme du code Wien2k 
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Chapitre 2 


Etude des propriétés structurales et 
électroniques des composés Ti1-:A:02 
(A= Ni, Ce et x =0, 0.25, 0.50, 0.75) 


Dans ce chapitre, nous avons utilisé l'approximation GGA pour explorer les propriétés 
structurales, électroniques de l’oxyde de titane TiO2 pur, dans sa phase rutile, et pour 


différentes concentrations des composés Ti1:.,Ax02 (A = Ni, Ce et x = 0.25, 0.50, 0.75). 
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L INTRODUCTION 


L’oxyde de titane TiO2 a suscité l'intérêt d'un grand nombre de chercheurs, en raison 
de ses excellentes propriétés, telles qu'une forte activité catalytique, une efficacité redox, une 
stabilité thermique et chimique élevée, il est non toxique... . Ces caractéristiques ont été 
largement utilisées dans plusieurs domaines de la recherche fondamentale et appliquée, citons 
par exemple le domaine des énergies renouvelables, la protection écologique de 


l'environnement, la photocatalyse..…. . . [1-4]. 


Dans ce chapitre, nous avons utilisé l'approximation GGA pour explorer les propriétés 
structurales et électroniques de l’oxyde de titane TiO pur, dans sa phase rutile, pris comme 
référence et pour différentes concentrations des composés T11.,Ax02 (A = Ni, Ce et x = 0.25, 


0.50, 0.75). 


IL DÉTAILS DE CALCUL 


Afin d'améliorer les propriétés optiques de l’oxyde de titane, dans sa phase rutile, les 
structures de bande, les densités d'état et les courbes d'absorption de TiO: pur et dopé au 
nickel et / ou au cérium ont été étudiées à l'aide du logiciel Wien2k [5]. Les différentes 
valeurs ont été calculées par la méthode des ondes planes augmentées et linéarisées à potentiel 
total (FP-LAPW) [6] dans le cadre de la théorie de la fonctionnelle de la densité DET [7]. Le 
potentiel de corrélation d'échange, de toutes les structures, est traité par l’approximation du 


gradient généralisée dans le paramétrage de Perdew-Burke-Erenterhof (PBE-GGA) [8]. 


TiO, rutile cristallise sous une structure tétragonale (Fig. 1) dans le groupe spatial 
P&;/mnm (N° 136 dans les tables internationales) Les paramètres de réseau sont : a = 4.593À 
et c = 2.959À [9]. Les états de valence pour les éléments utilisés sont les suivants: T1 (4s°3d?), 


O (2s°2p”), Ni (4s/3d°), et Ce (4f'5s°5p°5d'6s?). 


Pour étudier l'effet du dopage et l'effet du dopant A sur la géométrie structurale de 
TIO>, nous avons adopté une super cellule de 12 atomes dans laquelle la substitution de 
l’atome de titane Ti par A correspond à 25%, 50% et 75% de la concentration et il est 


remplacé par des atomes de Ni et Ce. 


Les rayons des sphères de muffin-tin (Rur) utilisés pour les atomes de Titane, 
Oxygène, Nickel et Cérium sont respectivement : 1.70, 1.60, 1.78 et 2.03 u.s. La valeur 1 
maximale des harmoniques sphériques étant de 8. Le nombre de fonctions de base FP-LAPW 


est choisi pour rendre R mr K max = 7. L'énergie de coupure qui définit la séparation du noyau 
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et des états de la valence est choisie égale à -6.0 Ry. L'espace réciproque est échantillonné par 
un maillage Monkhorst-Pack 8x8x6 avec suffisamment de 500 k-vecteurs dans la zone de 
Brillouin irréductible. Un critère auto-cohérent d'énergie totale avec une précision de 10* Ry 


est également utilisé. 


Figure 1 : Structure cristalline de l’oxyde de titane TiO, dans sa phase rutile. 
II. ETUDE DE L’OXYDE DE TITANE PUR 
IIL.1. Optimisation de la structure 


Les propriétés structurales de TiO; pur sont obtenues par une minimisation de l'énergie 
totale en fonction du volume de TiO; dans sa structure tétragonale (Fig. 2). Nous calculons les 
constantes de réseau, le module de compressibilité selon l'équation d'état de Murnaghan. 


BoV [(Vo/V)% BoVo 
B; | B-1 B,-1 


E(V) = Eo + (1) 


Où Fo et Vo sont l'énergie et le volume à l'équilibre. B et B' sont le module de compressibilité 
et sa dérivée. 

Les résultats de calcul obtenus, en utilisant l’approximation GGA, des constantes de 
réseau, du module de compressibilité, de sa dérivée ainsi que l'énergie totale sont résumés et 
comparés avec certaines données expérimentales et théoriques, disponibles de la littérature 


dans le tableau 1. 
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Figure 2 : Energies totales en fonction du volume pour TiO; rutile 


Tableau 1 Constantes de réseau a, c et module de compressibilité B de TiO; dans sa phase 
rutile. 


Nos calculs Autres calculs Expériences 
a (À) 4,634 4,630 [3]- 4,670[11] | 4,593 [10] 4,592 [12] 
C (À) 2,980 2,982 [3] -2,970[11] 2,958 [10] — 2,958 [12] 
B (GPa) 232,716 237,92 [3] 216 [10]-211 [12] 


Les constantes de réseau optimisées a= 4,634 À et c= 2,980 À pour TiO; pur sont en 
bon accord avec les données théoriques antérieures et les résultats expérimentaux tels que 


publiés dans les références suivantes [10, 12]. De même pour le module de compressibilité B. 


60 


Ceci montre le bon chois de la méthode de calculs. 


Le paramètre de maille a calculé est 1,8% plus grand que la valeur expérimentale, 
mettant en évidence que la GGA surestime la constante de réseau [13]. 


IIL2. Structure électronique 


La figure 3 montre les structures de bande électronique de TiO; en utilisant 
l'approximation GGA. Celle-ci montre que TiO: est un semi-conducteur à gap directe. 
L'énergie de la bande interdite calculée est de 1,81 eV au point l de symétrie élevé. Elle est 
inférieure à la valeur expérimentale [14]. La sous-estimation de la bande interdite est 
principalement décrite par le fait que le calcul par GGA ou le calcul par l'approximation de la 
densité locale (LDA) a la limitation de reproduire correctement les états électroniques 


inoccupés [15]. 


Energy(eV) 


Figure 3 : Structure de bande de TiO: rutile pur. 


Les densités d’états totales et partielles (TDOS et PDOS) de TiO, pur sont illustrées 
sur la Figure 4. La bande de valence (VB) est constituée essentiellement par l'état O 2p, alors 
que la bande de conduction (CB) est majoritairement constitué par l'état Ti 3d. Le passage des 
électrons des états O 2p aux états Ti 3d contrôle le début de la réaction dans l'absorption des 


ultraviolets. 
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Figure 4: TDOS (a) et PDOS (b) de TiO, pur dans sa phase rutile. 
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IV. ETUDE DES PROPRIETES STRUCTURALES ET ELECTRONIQUES DES 
COMPOSES Ti:.xA;02 (A= Ni, Ce et x = 0.25, 0.50, 0.75) 


IV.1. Introduction 


L’oxyde de titane est fortement dopé par des impuretés atomiques. C’est pourquoi, il 
existe des efforts de recherche académiques considérables pour contrôler le dopage de TiO; 


afin d’améliorer ses propriétés physiques existantes et prédire de nouvelles applications. 


Dans cette partie, nous avons aussi utilisé l'approximation GGA pour calculer les 
propriétés structurales et électroniques des composés Ti: ,A,02 (A= Ni, Ce et x = 0.25, 0.50, 
0.75) 


IV.2. Etude des propriétés structurales et électroniques des composés Ti.,CexO2 
(x = 0.25, 0.50, 0.75) 


IV.2. a. Propriétés structurales 


Dans ce cas, nous avons substitué 25%, 50% et 75% des atomes de titane T1 par des 
atomes de cérium Ce dans la matrice de l’oxyde de titane TiO, (Fig. 5). Nous avons effectué 
une optimisation par une minimisation de l'énergie totale en fonction du volume de TiO; dopé 


au Ce dans la structure tétragonale pour les différentes concentrations (Fig. 6) 


Figure 5 : Structure cristalline des composés Tii,CexO2 (x= 0, 0.25, 0.50, 0.75) 
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Figure 6: Variation de l’énergie totale en fonction du volume pour les différentes 


concentrations en Ce 
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Les résultats d'optimisations structurales pour TiO; dopé au Ce sont résumés dans le 
tableau 2. Les constantes de réseau (a et c) calculées, pour les composés Ti:,Ce;O», 
augmentent nettement avec la concentration du cérium lorsque l’oxyde de titane est dopé (Fig. 
7 et 8) alors que le module de compressibilité B diminue (Fig. 9).Ceci est certainement relié 


au rayon atomique du cérium qui est plus grand que celui du titane. 


Tableau 2 : Constantes de réseau (a et c) et modules de compressibilité B des composés 


TixCexO 
TiO: dopé Constantes de réseau Constantes de réseau modules de 
Ce % (a) (A) (c) (A) compressibilité B 
25% 4,796 3,199 185,835 
50% 4,850 3,495 161,690 
75% 4,901 3,792 147,129 


IV.2. b. Propriétés électroniques 


Les calculs des propriétés électroniques des composés Ti: ,Ce;O2 sont également 


effectués par l’approximation GGA. Les résultats de l’énergie de la bande interdite, pour les 
différentes concentrations (x = 0, 0,25, 0,50, 0,75) sont énumérés dans le tableau 3. Ceux ci 
montrent que les composés sont des semi-conducteurs et qu'ils présentent un écart de bande 


directe au point [de symétrie élevé. 


Tableau 3 : Les gaps énergétiques directs Eg(eV) des alliages Ti; ,Ce;O2 par rapport aux 


autres calculs théoriques et expérimentaux. 


Ti-xCexO2 Notre travail Exp Autres travaux 

0 1,81 3.00 [16] 1,56[3] 1,88[17] 
0,25 1,55 - - 
0,50 1,40 ; 
07 1,22 à : 
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Figure 7 : Variation du paramètre du réseau a deTi; ,Ce,;O; en fonction de la concentration x. 
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Figure 8 : Variation du paramètre du réseau c de Ti ;Ce;O en fonction de la concentration x. 
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Figure 9: Variation du module de compressibilité de Ti, ,Ce;O; en fonction de la 
concentration x. 
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La figure10 montre l’allure de la courbe reliant la bande interdite et la composition en 


cérium. On peut noter que la bande interdite varie avec la proportion d'impureté. On constate 
que Eg diminue avec la concentration en Ce. 
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Figure 12 : Densités d’états totales et partielles de Tio.75Ce25O, 
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Les structures de bandes électroniques, les densités d’état totales et partielles de 
Ti: xCexO par l’approximation GGA sont illustrées sur les figures 11, 12 et 13 
respectivement. La première figure montre que le maximum de la bande de valence (VBM) 
est presque sous le niveau de Fermi aux mêmes endroits que TiO, pur, tandis que le minimum 
de la bande de conduction (CBM) s’est déplacé vers le bas. En conséquence, la bande 


interdite du rutile TiO: dopé au Ce a diminuée. 
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Figure 1 3 : Densités d’états totales et partielles de Tio 50Ce500 
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Ceci est cohérent avec les résultats des densités d’états où la BC est constituée 
principalement par les orbitales Ce 4f et Ti 3d. Cette dernière est décalée vers les basses 
énergies par rapport à Ti1O2 pur suivant que la concentration en cérium augmente Ce 


phénomène est surtout l'effet de l'orbite Ce 4f à l'intérieur de la bande de conduction qui 


renvoi l'orbite Ti 3d vers le bas. 
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Figure 14 : Densités d’états totales et partielles de Tio 25Ceo 7502 
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IV.3. Etude des propriétés structurales et électroniques des composés Ti:.,Ni;O2 
(x = 0.25, 0.50, 0.75) 
IV.3.a. Propriétés structurales 


Dans ce cas, nous avons substitué 25%, 50% et 75% des atomes de titane T1 par des 
atomes de nickel Ni et effectué l'optimisation des paramètres structuraux. Ceux ci sont 
obtenus par une minimisation de l'énergie totale en fonction du volume de TiO: dopé au Ni 


dans la structure tétragonale (Fig 15). 


Les résultats d'optimisation sont résumés dans le tableau 4. Les constantes de réseau (a,c) et le 
module de compressibilité (B) calculés pour les composés Ti: ,Ni,0: diminuent nettement 
lorsque l’oxyde de titane est dopé par le nickel. Ceci est attribué à Ro, la différence de rayons 


covalents de Ni** en comparaison avec Ti**. 


Tableau 4: Constantes de réseau a, c et module de compressibilité B des composés Ti1.xNixO» 


TiO,; dopé au Ni % | Constantes de réseau | Constantes de réseau module de 
(a) (À) (c) (À) compressibilité B 
25% 4,585 2,959 223,053 
50% 4,521 2,924 220,073 
75% 4,466 2,888 213,689 


Les figures 16 et 17 représentent les paramètres de maille a et c des composés Ti1. 


xNixO2 pour différentes concentrations. Celles-ci montrent que lorsque la concentration en 
nickel augmente, les paramètres du réseau a et c diminuent. Le même phénomène est observé 
pour le module de compressibilité B calculé pour Ti:.,Ni;O2 qui diminue avec l'augmentation 


de la concentration x en Ni, comme le montre la Figure 18. 
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Figures 15 : Variation de l’énergie totale en fonction du volume pour les composésTi; ,Ni;O> 
(x = 0.25, 0.50, 0.75) 
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Figure 16 : Variation du paramètre du réseau (a) en fonction de la concentration x. 
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Figure 17 : Variation du paramètre du réseau (c) en fonction de la concentration x. 
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Figure 18 : Variation du module de compressibilité en fonction de la concentration x. 
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1V.3.b. Propriétés électroniques 


Les structures de bandes électroniques, les densités d’état totales et densités partielles 


de Tii.,NixO) sont montrées sur les figures 19, 20, 21 et 22 respectivement. 
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Figure 19 : Structures de bande des composés Ti: ,NixO. 


Nous observons l’existence de nouvelles bandes d’énergie (IELs : impurity energy 


levels) au dessus du niveau de Fermis. Celles-ci sont d’autant plus importantes que la 


concentration en nickel croit (Fig.19). 
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Figure 20 : Densités d’états totales et partielles de Tio 75Nio 250. 
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Ces résultats sont en bon accord aves les résultats des densités d’états qui montrent 
que les états d'impureté sont situés au-dessus du maximum de la bande de valence au dessus 
du niveau de Fermi et sont composés majoritairement par les états Ni 3d et O 2p. L’orbital Ni 


3d contribue à la bande de conduction en plus de l’orbital Ti 3d. 
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. Figure 21 : Densités d’états totales et partielles de Tio sNio.50O2. 
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Les électrons peuvent alors, être excités entre ces niveaux d'impuretés et également 
avec la bande hôte. Ainsi, les électrons Ni 3d contribuent de manière significative à la bande 
de valence. Ils peuvent être excités par les IELSs, ensuite, ils excitent le CBM par irradiation de 


la lumière visible 
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Figure 22 : Densités d’états totales et partielles de Tio25Nio 7502. 
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V. CONCLUSIONS 


En résumé de ce chapitre, nous avons étudié les propriétés structurales et électroniques 


de TiixAxO2 (A = Ni, Ce et x =0, 0.25, 0.50, 0.75) en utilisant l'approximation GGA 


Nos résultats concernant les propriétés structurales sont en bon accord avec les valeurs 
théoriques et expérimentales de la littérature. Nous avons étudié la dépendance des constantes 
de réseau, du module compressibilité et la bande interdite des composés Ti1.,Ax02 (A=Ni, Ce) 
en fonction de la concentration en Ce et en Ni. Nous avons observé l'effet de (A=Ni, Ce) sur 
les propriétés électroniques de Ti, ,A;O> pour (x = 0, 0,25, 0,50, 0,75).Les structures de 
bandes calculées et les densités d'états indiquent un caractère semi-conducteur des matériaux 


impliqués pour Ti:.,Ce;O», et un caractère métallique des matériaux impliqués de Ti:.,NixO. 
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Chapitre 3 


Etude des propriétés structurales 
électroniques et optiques du composé 
Ti0.50Ce0.25N10.2502 


Dans ce chapitre, nous avons utilisé l'approximation GGA pour calcules les propriétés 
structurales, électroniques et optiques du composé Tio50Ce025Nio2502de l’oxyde de titane 
avec TiO;, dans sa phase rutile, pris comme référence. 


Chapitre 3 : Etude des propriétés structurales, électroniques et optiques des 
composés Tis 59C€0 25 Ni0. 2502 
L. INTRODUCTION 
II. PROPRIETES STRUCTURALES 
III. PROPRIETES ELECTRONIQUES 
IV. PROPRIETES OPTIQUES 
IV. 1. Introduction 
IV.2. Fonction diélectrique 
IV.3. Indice de réfraction et coefficient d'extinction 
IV.4. Coefficient d'absorption 


V. PREDICTION DU POTENTIELS REDOX DE L’ELECTRONEGATIVITE 
ATOMIQUE 
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L Introduction 


Dans ce chapitre, nous avons utilisé l'approximation GGA pour calculer la structure, 
électronique et l'absorption optique de l’oxyde de titane TiO2 pur,dans sa phase rutile, pris 
comme référence et co-dopé au nickel (Ni) et au cérium (Ce). L'analyse théorique fournie 
peut être une explication possible du phénomène de l'absorption optique pour le composé 


(Ni,Ce)TiO». 
IL. Propriétés structurales 


Les propriétés structurales des composés Ti05, Ti3N1Og, T13CeO% et TiNiCeO;(Fig. 
1)sont obtenus par une minimisation de l'énergie totale en fonction du volume de TiO:et de 
TiO; dopé au nickel et / ou au cérium dans la structure tétragonale (Fig. 2).Nous calculons les 
constantes de réseau, le module de compressibilité en ajustant l'énergie totale en fonction du 


volume selon l'équation d'état de Murnaghan. 


E(V) = Eo + 


Bo 
Fe +1] BoVo «) 


Bel Bi _Bi-1 


Où Eo et Vo sont l'énergie et le volume à l'équilibre. B et B'sont le module de compressibilité 


et de sa dérivé. 


Les propriétés structurales à l'équilibre telles que les constantes de réseau et le module 
de compressibilité sont calculées en utilisant l’approximation GGA. Les résultats obtenus sont 


résumés dans le tableau 1. 


Tableau 1 Constantes de réseau a, c, module de compressibilité B des composées Ti3N1Os, 


Ti3CeOs et Ti NiCeOs 


Structure (a) (À) (c) (À) B (GPa) 
TiO> 4,634 - 4,593[1] 2,980 2,958 [1] 232,716 216 [1] 

4,630[2] 2,982 [2] 237,92[2] 

Tio35Nio 2502 4,585 2.959 223,033 

Tio75Ce02502 4,796 3,199 185,835 

Tio.sNio25Ce0 2502 4,772 3,171 171,448 
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(c) (d) 


Figurel : Structures cristallines de (a) TiO2pur, (b) TiO2 dopé Ni (Tio75Nio25 On), (c) TiO2 
dopé Ce (Ti 75 Ceo.25 O2),(d) TiO2co-dopé Ni-Ce (Ti0,50 Nio,25 Ceo 250). 

En comparaison avec les paramètres structuraux purs de TiO», la constante de réseau 
cristallin (ao, Co) diminue légèrement (1% et 0,7% respectivement) avec le TiO,dopé au Ni 
pour la structure Ti;N104. Cependant il augmente (de 3% et 7% respectivement) avec le rutile 
Ti dopé au Ce pour la structure Ti3CeOs. Pour le rutile TiO2co-dopé (Ni, Ce), une 
augmentation de 2% et 6% respectivement a été trouvée pour la structure Ti: NiCeOs.Ceci est 
attribué à Ro, la différence de rayons covalents de Ni**, Ce** en comparaison avec Ti* (R Ni* 


#2 0,48 À, R Ce**= 0,97 À, R Ti‘* = 0,605 À). 


IL. Propriétés électroniques 


Les résultats de calculs des structures de bandesdes composés TiO2,T13NiOg, Ti3CeOs 
et Ti2NiCeOs sont représentés sur la figure 3. Il est bien clair sur la figure 3a que le maximum 
de la bande de valence et le minimum de la bande de conduction présentent un écart 
énergétique direct au point ['. L'estimation de cette bande interdite pour TiO: pur est de 1.81 
eV. Cette valeur est inférieure à la valeur mesurée expérimentalement, en raison de la 
discontinuité du potentiel de l'échange et de corrélation, non pris en compte dans le cadre de 


la DET. 
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Figures 2 : Variation de l’énergie totale en fonction du volume pour 
a) Ti0,75C€0,25 O2 b)Ti0,75Ni0,2502 C)Ti0,50 Nio,25 C€o 2502. 
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Cependant, la bande interdite indirecte calculée pour TiO, rutile dopé au Ni est 
d'environ 2,25 eV, comme indiqué sur la figure 3b. Celle ci montre l'existence de deux 
niveaux d'énergie d'impureté (IELs). Ils sont situés au-dessus du maximum de la bande de 
valence. Ce résultat indique que des électrons peuvent être excités entre ces niveaux 
d'impuretés et également vers la bande hôte. La figure 3c montre une bande interdite indirecte 
de 1,55 eV pour le rutile dopé au Ce, qui est inférieure à la valeur du TiO pur. Finalement 
pour TiO; co-dopé (Ni, Ce) (Fig.3d), le gap calculé est d'environ 2,03 eV. Le co-dopage du 
nickel et du cérium peut ajuster la structure de bande du dioxyde de titane rutile, en plus il a 


développé des niveaux d'énergie d'impureté dans la bande interdite. 
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Figure 4 : Densités d’états totales et partielles de TiO:pur. 
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Les figures4, 5, 6 et 7 présentent les courbes des densité d'états totales et partielles 
calculée pour tous les composés. Pour TiO: pur, la bande de valence (VB) est 
fondamentalement constituée par l'état O 2p, bien que la bande de conduction (CB) soit en 
grande partie formée par l'état Ti 3d, comme montré sur la figure.4b. Cela suppose que le 
passage des électrons des états O 2p à Ti 3d est responsable de l'apparition de l'absorption des 
ultraviolets. 
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Figure 5 : Densités d’états totales et partielles des composés Ti3NiO». 
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Pour l'oxyde de titane dopé au Ni (Fig.S), les états d'impureté sont situés au-dessus du 
maximum de la bande de valence, composés des états O 2p et Ni 3d. Il est évident que les 
électrons Ni 3d ont une contribution significative à la bande de valence. Les électrons du 
VBM peuvent être excités aux IELs puis excités aux CBM par irradiations de la lumière 
visible. Pour le rutile dopé au Ce (Fig.6), la courbe est semblable à celle de TiO: avec 
diminution du gap. Pour TiO; co-dopé (Ni-Ce) (Fig.7), la bande interdite directe a augmenté 
par rapport à Ti10, mais elle a diminué par rapport à TiO dopé au Ni.Ce phénomène est 


notamment l'effet de l'orbite Ce 4f à l'intérieur de la bande de conduction. 
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Figure 6 : Densités d’états totales et partielles des composés Ti3CeO:. 
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Ce phénomène a été provoqué par la bande de valence, qui est principalement 
constituée par les états Ni 3d et O 2p et la bande de conduction dominée par les états Ce 4f 
localisé au milieu de celle ci. De plus, les niveaux d'énergie d'impureté développés par les 
états Ni 3d sont situés en haut de la VB. Dans ce cas, les électrons dans la bande de valence 
peuvent être excités à la bande de conduction après avoir traversé les IELS par l'absorption de 
la lumière visible. Les IELSs peuvent être considérées comme des centres de pièges pour des 
paires d'électrons trous photo-excités, ce qui peut également réduire le taux de recombinaison 
des porteurs de charge et prolonger la durée de vie des porteurs. Par conséquent, ces niveaux 
d'énergie d'impureté ont des avantages dans l'extension de la longueur d'onde sensible à la 


région de la lumière visible. 
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Figure 7 : Densités d’états totales et partielles des composés Ti,NiCeO; 
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IV. Propriétés optiques 
IV.1. Introduction 


De manière générale, les propriétés optiques sont décrites à l'aide de la fonction 
diélectrique complexe £& (w) qui est donnée par & (©) — £1 (@) + 162 (œ). Les autres paramètres 
optiques tels que l'indice de réfraction n(o), le coefficient d'extinction k(w), le coefficient 


d'absorption, la réflectivité et la transmittance peuvent être déduits directement de £(œ) 
IV.2. Fonction diélectrique 


La fonction diélectrique £ (@) — £1 (@) + 1£2 (©). où e1 (©) est la partie réelle 
quireprésente la dispersion des photons incidents par les matériaux tandis que £> (@) est la 
partie imaginaire qui résulte des transitions inter-bande entre les états occupés au-dessous du 


niveau de Fermi et les états inoccupés au niveau de Fermi. Ceci décritl'absorption de photons. 


La partie imaginaire diélectrique est donnée mathématiquement comme suit 
4T?e? nu. 
Ex(@)ag = nr D feririelrslr) Wi(1 —W,)ô(E; — E; — ho)d#k (2) 
if 


où(f|P,li}et (ilPe| f }sont les éléments de matrices dipolaire, f et i sont les états finals et 


initiales respectivement.,; (resp. W£ ) est la fonction de distributions de Fermi centrée en 


E;(resp. Er ). 


La partie réelle £, (w) est calculée à partir de &> (©) en utilisant la relation de Kramers- 


Kronig [3] donnée par : 


2 F &'£(&') 
re 


dw' (3) 


où P et « sont la constante d'intégration et la fréquence des photons incidents. 


La figure 8 montre £, (@) et #2 (©) de TiO, pur et du TiO; dopé au nickel (Ni) et / ou au 


cérium (Ce). 
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TiO, pur 

——— Ni doped TiO, 
Ce doped TiO, 
Ni-Ce co-doped TiO, 


Dielectric real part 


Energy (eV) 


——— TO, pur 

——— Ni doped TiO, 
Ce doped TiO, 
Ni-Ce co-doped TiO, 


dielectric imaginary part 


Energy (eV) 


Figure 8 :Fonctions diélectriques de TiO; pur et du TiO; dopé au nickel (Ni) et / ou au 
cérium (Ce): (a) partie réelle et (b) partie imaginaire. 
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Nous observons (Fig. 8a) une augmentation significative, de la partie réelle £, (w)de la 
fonction diélectrique, après dopage avec (Ni) et / ou au cérium (Ce), aux faibles énergies. La 
fonction diélectrique statique, £1(0) correspondant à une pulsation nulle, pour TiO:, TiO2 dopé 
au nickel (Ni), TiO; dopé au cérium (Ce) et TiO,co-dopé au nickel et au cérium (Ni + Ce) 


sont respectivement 7,40, 24,14, 7,80 et 19,25. 


La fonction diélectrique imaginaire £2 (@) des différents composés dans l'intervalle 0-6 
eV est illustrée sur la figure 8b. La courbe de TiO: rutile montre clairement un bord à environ 
1,87 eV, qui provient principalementdes déplacements des électrons entre les états O 2p de 
VBM et Ti 3d du CBM. En générale, cette valeur d'énergie est appelée gap optique. Les pics à 
2,80, 2,18, 2,76 eV sont principalement dues aux déplacement entre les états d’impureté 
occupés peu profonds et les états vacants au minimum de la bande de conduction, pour Ti0: 
dopé Ce, dopé Ni et co-dopé (Ni+Ce). Ces énergies sont dans le champ des photons visibles. 
Ces valeurs pourraient être liées à une bande interdite décalée dans le système TiO, rutile 
dopé. Ce déplacement de la constante diélectrique vers les énergies plus faible contribue à 


l'amélioration de l'absorption de la lumière visible. 
IV.3. Indice de réfraction et coefficient d'extinction 


L'indice de réfraction est un paramètre physique très important lié aux interactions atomiques 


microscopiques Il peut être décrit par: 


l£2(w)| + (wo) 


, @) 


n(w) = 
Le coefficient d'extinction définit comme une mesure du taux de diminution de la lumière 
transmise par diffusion et absorption pour un milieu. Il est écrit sous la forme suivante: 


e}(w) 


ROSES 2e (5) 


Où #1 (@) et £2 (w) sont la partie réelle et la partie imaginaire diélectrique. 


La relation entre la fonction diélectrique et l'indice de réfraction complexe N = n + ik est 


donné par les relations : £1(@) = n°- k? où n = &> (m}/k. 


En particulier, la valeur de n pour la fréquence nulle no est donné par : 
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mn Li0 pur 


——— Ni doped TiO, 
Ce doped TiO, 
Ni-Ce co-doped TiO, 


Refractive Index 


——— Ni doped TiO, 
Ce doped TiO, 
——— Ni-Ce co-doped TiO, 


Extinction Coefficient 


Energy (eV) 


Figure 9 : (a) Indice de réfraction et (b) Coefficient d'extinction, de TiO: pur et de TiO: dopé 
au nickel (Ni) et / ou au cérium (Ce). 
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ni [ Oe (6) 


&w 
La Figure 9. Illustre l'indice de réfraction calculé n (&) et le coefficient d'extinction k 
(wo). Les indices de réfraction statiques déduits à partir des courbes de la figure 9 sont 2,70, 
5,15, 2,79 et 4.39 pour TiO: pur, TiO2 dopé (Ni), dopé(Ce) et co-dopé(Ni + Ce) 


respectivement. 
IV.3. Le coefficient d'absorption 


Le coefficient d'absorption est une mesure du taux de diminution de l'intensité du 
rayonnement électromagnétique (en tant que lumière) lorsqu'il passe à travers une substance 


donnée. Il est donné symboliquement comme : 


1/2 
a(w) = V2 w | [e2(c) + E2(w) — ex) (7) 


Les spectres d'absorption de TiO: pur et dopé sont montrés sur la figure 10. Le seuil 
d'absorption de TiO, rutile pur est limité à la zone de lumière ultraviolette seulement. 
Néanmoins, pour le TiO; rutile dopé, ces bords d'absorption se déplacent vers la lumière 
visible après dopage au nickel ou au cérium. Le phénomène ci-dessus est encore plus évident 
dans le cas de TiO; rutile co-dopé (Ni, Ce), qui montre une absorption de la lumière visible 
plus intense et plus élevée que TiO: pur ou simplement dopé. Ces résultats permettent de 
mettre en évidence que le co-dopage de dioxyde de titane avec Ni et Ce améliore l’absorption 


de la lumière visible. 
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TiO2 pur 

Ni doped TiO2 

Ce doped TiO2 
Ni-Ce co-doped TiO2 


Absorption (a.u.) 


300 400 500 600 700 800 
(nm) 
Figure 10 : Coefficient d'absorption du TiO,; pur et du TiO, dopé au nickel (Ni) et / ou au 


cérium (Ce) 
V. Prédiction des potentiels redox de l’électronégativité atomique 


L'électronégativité est l'un des paramètres atomiques les plus fondamentaux qui 
exprime numériquement la tendance à attirer les électrons vers les atomes dans les cristaux. 
Parr et al. avait donné la définition précise à l'électronégativité d'un atome neutre selon la 


DFT [4], qui est une valeur négative du potentiel chimique dans l'état fondamental: 
: : ( ) à 
ATH x) (8) 


Où E et N étant l'énergie totale électronique de l'état fondamental et le nombre d'électrons, 
respectivement. y et 1 étant l'électronégativité et le potentiel chimique de l'atome. Aïnsi, la 


formule (8) peut être exprimée par une équation différentielle, correspondant à l’équation 


suivante : 
ÔE 
= (5) = [E(N + 1) — E(N) + ECN) — E(N — 1)]/2 & (1 + A)/2 (9) 


Dans cette équation, I est l'énergie d'ionisation et A est l'affinité électronique pour l'atome 


[SI.II se trouve que la formule (9) a la même expression avec la définition de Mulliken de 
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l'électronégativité [6]. Selon la formule (9), nous avons calculé les électronégativités des 


atomes O, Ti, Ce et Ni.lesvaleurs sont résumées dans le tableau 2. 


Tableau 2 : Energie d'ionisation (1), affinité électronique (A) et électronégativité (x) des 


atomes. 
Atomes I(eV) [7] A(eV) [7] x(eV) 
O 13.618 1.451 7.540 
Ti 6.820 0.079 3.450 
Ce 5.538 0.650 3.094 
Ni 7.635 1.156 4.396 


L'électronégativité de Mulliken d'un composé semi-conducteur peut être définie 
comme la moyenne géométrique des électronégativités des atomes constitutifs [8,91, 


correspondants à la formule suivante : 


1 
= CT (10) 


Dans cette formule y, n et N représentent l'électronégativité de l'atome constituant, le nombre 
d'espèces atomiques et le nombre total d'atomes dans le composé semi-conducteur, 


respectivement. 


Nous avons calculé les électronégativités de TiO: pur, de TiO: dopé au Ce, de TiO2 dopé au 
Ni et de TiO: co-dopé (Ce + Ni) sur la base de la Formule (10), comme indiqué dans le 


Tableau 3. 


L’électronégativité de Ti0: dopé au Ce est légèrement inférieure à celle deTiO: pur, 
tandis que celles deTiO; dopé au nickel et co-dopé (Ce + Ni) sont légèrement supérieures. 


Ceci est en relation avec l'électronégativitéde la nature de chaque élément où :4YNi>Ti>YCe. 


Les résultats calculés montrent que les trois composés dopés et co-dopé possèdent une 
plus grande capacité réductrice liée à l'augmentation deleurs électronégativités.En outre ceci 
prouve que le système TiO2co-dopé (Ce + Ni) a la capacité de réduire H*pour produire de 


l'hydrogène H. 
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Tableau 3 : Electronégativité x(s)des systèmes TiO) pur, de TiO; dopé au Ce, de TiO; dopé 
au Ni et de TiO,co-dopé (Ce + Ni). 


Systèmes TiO;pur TiO,; dopé Ce TiO,; dopé Ni TiO;co-dopé (Ce+Ni) 


x(eV) 5.810 5.763 5.928 5.880 


De plus, suivant les calculs des structures électroniques et des propriétés optiques, 
nous savons que les états hybrides Ce 4f et Ni 3d peuvent réduire la recombinaison des paires 
d'électron-trous photo-générés dans le système TiO,co-dopé(Ce + Ni) et améliorer 
l'absorption optique par déplacement du seuil d'absorption de la lumière UV à la lumière 


visible par rapport au système TiO: pur. 


La figure 11 montre le mécanisme physique et chimique de la dissociation 


photocatalytique de l'eau pour générer de l'hydrogène dans le système T10, (Ce + Ni). 


Par conséquent, le système TiO: co-dopé (Ce + Ni) pris comme photocatalyseur peut 


présenter une excellente activité photocatalytique dans le domaine de la lumière visible pour 


la production d'hydrogène par dissociation de l'eau par l'effet synergique du co-doping (Ce + 


Dissociation photocatalytique de l'eau 


Ni). 


Lumière 
du soleil 


Génération d'hydrogène 


Figure 11 : Illustration schématique du mécanisme physique et chimique de la dissociation 
photocatalytique de l'eau pour la production d'hydrogène dans le système T10: (Ce + Ni). 
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CONCLUSION 
En résumé de ce chapitre, une étude théorique de la structure, électronique et les propriétés 


optiques du TiO; pur et du Ti0: dopé au nickel (Ni) et / ou au cérium (Ce) ont été rapportées. 
Cette étude a été réalisée à l'aide de DET avec approximation GGA. On trouve que la bande 
interdite d'énergie calculée et le paramètre de réseau de TiO; pur sont proches des paramètres 
expérimentaux et en bon accord avec d'autres calculs théoriques. Il est également montré que 
l'incorporation de nickel (Ni) et / ou au cérium (Ce) dans TiOzaffecte considérablement les 
propriétés électroniques et optiques par rapport au TiO> pur.Par exemple, l'écart d'énergie 
optique a augmenté de 12% sous le co-dopage Ni+Ce. De la fonction diélectrique imaginaire, 
nous avons établi que le décalage vers le rouge est également réalisé sous le dopage cérium 


(Ce) indiquant l'importance de ses états occupés par le donneur 4f. 
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Conclusion générale 
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Ce travail de thèse constitue une étude des propriétés structurales, électroniques et 
optiques de l’oxyde de titane TiO: dans sa phase rutile dopé par le nickel ou le cérium. Cette 
étude est basée sur une méthode ab-initio dite méthode des ondes planes linéairement 
augmentées à potentiel total (FP-LAPW) implémentée dans le code Wien2k et dans le cadre 


de la théorie de la fonctionnelle de la densité, en utilisant l’approximation GGA. 


Les constantes de réseau calculées pour TiO: pur sont très proches des résultats 
expérimentaux, ce qui prouve que nos paramètres de calcul sont valides. La valeur calculée 
est de 1% supérieur à la valeur expérimentale, mettant en évidence la surestimation par la 
GGA de la constante de réseau. De même, pour les autres grandeurs calculées de TiO; pur, 


sont en accord avec les études théoriques et expérimentales de la littérature. 


Nous avons également effectué des calculs sur la structure électronique, en 
déterminant les structures de bandes électroniques et les densités d’états électroniques totales 
(TDOS) et partielles (PDOS) qui ont montré que Tio75N102502 est un conducteur et les deux 


composés, T10: et Tio75Ce02502 sont des semi-conducteurs. 


Une étude sur le composé Ti1o,50Ce025Ni02502 a été présentée dans la deuxième partie. 
Nous avons calculé le paramètre de réseau, le module de compressibilité, les gaps 


énergétiques et les grandeurs optiques. 


A la fin, les résultats montrent que le Ti02 dopé au Ni induit des niveaux d'énergie 
d'impureté supérieurs à la bande de valence, qui peuvent être responsables de l'absorption de 
la lumière visible, tandis que le dopage avec Ce rétrécit la bande interdite du T1: rutile. La 
combinaison de Ni et de Ce comme co-dopant de TiO; améliore significativement l'activité 


photocatalytique du TiO; rutile dans la gamme de la lumière visible. 


101 


